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Abstract 


The contents of this report covers: (i) development of optimal geometry 

for crowned spur gears; (ii) methods for their generation; and (iii) tooth 
contact analysis (TCA) computer programs for the analysis of meshing and 
bearing contact of the crowned spur gears. 

The developed method for synthesis is used for the determination of the 
optimal geometry for crowned pinion surface and is directed to reduce the 
sensitivity of the gears to misalignment, localize the bearing contact, and 
guarantee the favorable shape and low level of the transmission errors. 

A new method for the generation of the crowned pinion surface has been 
proposed. This method is based on application of the tool with a surface of 
revolution that slightly deviates from a regular cone surface. The tool can 
be used as a grinding wheel or as a shaver. The crowned pinion surface can be 
also generated by a generating plane whose motion is provided by an automatic 
grinding machine controlled by a computer. 

The TCA program simulates the meshing and bearing contact of the 
misaligned gears. The transmission errors are also determined. 
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Summary 


Spur gears are widely used to transmit mechanical power between 
parallel axes of drive shafts. Namely, such gears are used in 

planetary trains of helicopters. It is well known that spur involute 
gears are very sensitive to gear misalignments. Misalignment will 
cause: the shift of the bearing contact toward the edge of the gear 

tootn surfaces and transmission errors that increase gear noise. The 
previous methods for crowning of pinion tooth surface have been 
directed at providing of a favorable bearing contact. The 

transmission errors of the misaligned gears have been ignored. The 
new approach developed oy the authors of this report is directed at 
the synthesis of spur involute gears with a crowned pinion surface 
tnat provides: (i) a parabolic type of low transmission errors for 

the misaligned gears and ( i i ) a localized bearing contact. New 

methods for the generation of crowned pinion surface proposed by the 
authors can be readily applied by the gear industry. 
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1. Introduction 


Spur involute years are very sensitive to gear misalignments that 
cause: (i) the shift of the bearing contact to the edge of the gear 

tooth surfaces and (ii) transmission errors that increase gear 
noise. Many efforts have been done to improve the bearing contact of 
i misaligned spur gears by crowning of pinion tooth surface. Wildhaber 
> ( ly62) has proposed various methods of crowning that can be achieved 

in the process of gear generation. Maag engineers (see References) 
have used crowning as (i) longitudinal corrections (Fig. 1.1a), (ii) 
modified involute tooth profile uniform across facewidth (Fig. 1.1b), 

(iii) combination of longitudinal correction and uniform modified 

i 

1 profile (Fig. 1.1c) and (iv) topological modification (Fig. l.ld) that 
can provide any deviation of the crowned tooth surface from a regular 
! involute surface. 

I 

( The main purpose of the previously proposed methods for crowning 

i 

I 

; was to improve the bearing contact of the misaligned gears. This 
i resulted in only part of the solution. The transmission errors of the 

misaligned spur gears (a main source of the noise), were ignored. The 

I 

influence of gear misalignment on the transmission errors has not been 
; investigated. Maag's method of topological modification does not 
i quantify the relationship between the surface deviation and the 

I 

I transmission error. Also the optimal geometry for the crowned pinion 
i surface has not been proposed. 

The contents of this report cover the solutions to the following 
problems: (i) optimal yeometry of a crowned pinion tooth surface; 

(ii) new method of crowning based on application of a tool with a 
surface of revolution (the tool surface is slightly deviated from a 
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cone surface) (iii) the development of Tooth Contact Analysis (TCA) 
programs for the determination of transmission errors for misaligned 
gears and their bearing contact. 

The development of the optimal geometry of a crowned pinion gear 
tooth surface is based on the following considerations: 

(i) Misaligned spur gears with a crowned ninion tooth surface 
can provide transmission errors A 2 ^ ^ i ^ °f two t VP es that are shown in 
Fig. 1.2a and 1.2b, respectively. The transmission errors are 
determined with the equation as 


A<f> 

2 ( $ i ) ^2^1^ 

N 1 

" N 2 

♦1 



(1. 

1) 

Here : 

and N 2 are 

the 

numbers of gear teeth; 

♦l 

and *2 

are 

the 

angles 

of gear rotation; 

4) 2 ( * ^ ) the function 

that relates 

the 

angles 

of rotation of 

gear 

and pinion if the pinion 

is 

crowned 

and 

the 


. . N 1 

gears are misaligned; 2 q = yj - <f>^ is the theoretical relation between 
the angles of rotation of the gears in the ideal case where the gears 
are not crowned, not misaligned and the transmission errors do not 
exist. Type 1 of transmission errors are not acceptable because the 
change of tooth meshing is accompanied with an interruption or 
interference of tooth surfaces. Type 2 of transmission errors is more 
preferable if the level of the transmission errors does not exceed a 
prescribed limit. 

(ii) On the first glance crowning should be directed at 
providing an exact involute shape in the middle cross-section 
(Fig. 1.3). In reality this type of crowning is not acceptable 
because the misaligned gears will transform rotation with transmission 
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errors of type 1 (Fig. 1.2a). This is the reason why the authors have 
decided to synthesize a specific crowned pinion tooth surface. Such a 
pinion can provide transformation of rotation with a parabolic type of 

transmission error function when it is in meshing with a gear with a 

i 

i regular involute tooth surface. This type of function of transmission 
errors is synthesized for ideal gears that do not have any 

j misalignment. Then, the tendency to provide parabolic transmission 
errors can be extended for the misaligned gears and the discontinuance 
of mesning can be avoided. It is necessary to emphasize that the 

proposed method of synthesis provides a shape in the middle cross- 
' section of the pinion tooth that deviates in a certain way from the 
i involute curve that is shown in Fig. 1.3. The longitudinal deviation 

i 

; from a straight line is not related with the transmission errors but 
I with the desired dimensions of the instantaneous contact ellipse for 
the gear tooth surfaces. The proposed pinion tooth surface can be 
generated by a plane as the chosen generating surface. The motions of 
j the plane with the respect to pinion must be controlled by a computer. 

i 

(iii) Another method of pinion crowning is based on application 
of a surface of revolution that slightly deviates from a regular tool 
j conical surface (Fig. 1.4). Such a tool can be used as a grinding 
wneel or as a shaver. The motions of the tool and the gear being 
generated are related similarly to the motions of a rack-cutter and 
the gear (see Chapter 4). A tool with a regular conical surface can 
generate a crowned pinion tooth surface whose middle cross-section 
represents an involute curve (Fig. 1.3). However, this type of a 
crowned pinion tooth surface is not desirable because the misaligned 
gears provide transmission errors of type 1 (Fig. 1.2a). This is the 
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reason why a tool with the surface of revolution is to be used instead 
of a conical surface. 

(iv) The evaluation of the bearing contact and transmission 

errors for the misaligned gears as well as the investigation of the 
influence of errors of gear assembly requires application of TCA 
programs. Such programs have been developed by the authors. The 

programs are based on the following algorithms: 

(a) The contacting gear tooth surfaces are represented in a 
fixed coordinate system, Sf, that is rigidly connected to the gear 
nousing (Fig. 1.5). 

(b) Tne continuous tangency of gear tooth surfaces is provided 

if tne position-vectors and surface unit normals for the contacting 
surfaces coincide at the contact point at any instant. Then, we are 
able to determine the path of contact on the gear tooth surfaces and 
tne relations between the angles of rotation of the output and input 
gears. Knowing function we can determine the deviations of 

this function from the prescribed linear function i.e. the 
transmission errors. 

(c) Due to the elasticity of the gear tooth surfaces the surface 

contact is spread over an elliptical area. The dimensions and 

orientation of the instantaneous contact ellipse depend on the 
principal curvatures and principal directions of the contacting tooth 
surfaces. The bearing contact is determined by the developed TCA 
program as the set of the contact ellipses that move over the 

contacting surfaces in the process of motion. 
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2. Optimal Pinion Geometry 
2.1 Basic Consideration 

( It is assumed that only the pinion will be provided with a 

l 

| crowned tooth surface while the gear is provided with a regular 
I involute surface. The crowned pinion tooth surface deviates from a 
regular involute surface. The topology of the crowned pinion tooth 
| surface must satisfy the following requirements: (i) the pinion tooth 

surface is a continuous one; (ii) the pinion and gear tooth surfaces 

t . 

| are in contact at a point at any instant; (in) the principal 

i directions and curvatures of the contacting surfaces must be related 

j 

i in such a way that the instantaneous contact ellipse will be of 
appropriate dimensions; (iv) the function of transmission errors of 
tne misaligned gears should be of a parabolic type and its level must 
not extend the prescribed limits. The sat isf ication to the described 

I 

I conditions is considered in the following sections of this part. 

i 

| 

! 2.2 Parabolic Type of Function of Kinematic Errors 

I Initially we will consider that the gears are not misaligned, the 

pinion tooth surface is crowned and the gear tooth surface is a 
[regular involute surface. The tooth surfaces are in point contact at 
I every instant and the points of contact lie in the middle cross- 
section. Thus, initially we consider tne meshing of gear tooth 
i surfaces as a planar gearing. However, when simulating the meshing of 
gears by TCA program, we will consider the mesh of them as spatial 
gears . 

Figure 2.2.1 snows the shapes of the gears which are in contact 
at a point that lies on the center distance. The shape of the gear 
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tooth profile is a regular involute curve, and the shape of the pinion 
tooth profile is to be determined. It is given that the shapes must 
provide transformation of rotation represented by the function (Fig. 
2.2.2a) 


4 2 ^ 1 



+ A<t>2 (<(>]_ ) 


( 2 . 2 . 1 ) 


where A 4> 2 ( ^ ) is a function of a parabolic type that represents the 

function of transmission errors; and N 2 are the numbers of gear 

teeth. We have to emphasize that the proposed method for synthesis 

provides low transmission errors for aligned gears and the function of 

transmission errors is a parabolic one. Then, the misaligned gears 

will keep the tendency to transform rotation with the same type of 

transmission errors and the discontinuance of contact by meshing will 

be avoided. This means that we have to avoid the appearance of 

transmission errors that are shown in Fig. 1.2a. Function 

N 1 

^ 2 ( ^ l ) cont ains a linear function represented by that relates 

the angles of rotation of the "ideal" input and output gears. Figure 
2.2.2a shows the function ^ 2 ( ^ 1 ) * The function of transmission 
errors is represented to a larger scale in Fig. 2.2.2b as a periodic 
function of a parabolic type. The derivative 


Atf) 2 U]^) 


d 

d <j)^ 


( A <f> 2 ( 4^) ) 


( 2 . 2 . 2 ) 


is shown in Fig. 2.2.2c. 


For a 


"pure" parabolic function this 
It is evident that the second 


derivative is a linear function 
derivative represented by 
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I 


d 2 

2 (A* 2 (* 1 )) 



( <^2 ( <j>i ) ) 


is negative. The transformation of rotation by the gears will be 

i 

performed with a non-constant gear ratio 


m 21 


“2 

0)1 



A<j>2 ( <f' 1 ) 


(2.2.3) 


[The instantaneous center of rotation I is moved along the gear center 
(distance in the process of gear rotation. Figure 2.2.1 shows the 

I 

instantaneous location of I. 

i 

I More details about the geometry of the pinion that corresponds to 

I 

the discussed function 4*0 < be given in Chapter 3. 

! 

I 


2 . 3 Relation Between the Principal Curvatures and Directions 

The relation between the principal curvatures and directions 

the contacting surfaces and the derivative of the gear ratio 

i 

Represented as follows (Litvin, 1987) 


of 

is 


b ll 

b 12 

b 13 


(2.3.1) 

b 21 

b 22 

b 23 

= 0 

b 31 

b 32 

b 33 




For the case of spur gears with the crowned pinion tooth surface the 
slements of the determinant of Eq. (2.3.1) are represented by the 
rollowing equations 
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(2.3.2) 


b 21 = D 12 


b ll b 12 


L- b 12 b 22 


'31 


= b 


- K-, 


13 

( 1 ) 


+ K 


32 

( 2 ) 


= b 


23 


- K 


( 1 ) 

II 


+ 1C 


( 2 ) 

II -i 


(2.3.3) 


' b 13 " 


r * (i 2 ) n 

-e T • a) 

rv /v 

- b 23 - 


- (12) 

L J 




(2.3.5) 


+ ( ) 2 m' (nxk 0 ) . ( rl 1 ^ -C ) 

z i ~ ~z ~r ~ 


Figure 2.3.1 shows the tangent plane to the gear tooth surfaces at the 
point of contact I (see Fig. 2.2.1). The principal directions are 
represented by the unit vectors e T and e TT that are the same for the 
Doth contacting surfaces. Vector e^ is shown in Fig. 2.3.1 and is 
represented in coordinate system Sf as follows 


sin^ c 

-cos lb 

v c 

0 


(2.3.6) 


wnere tp is tne pressure anyle at point I. Vector e TT is directed 
C 1 

parallel to the negative Zf-axis and is represented by 


3 


(2.3.7) 



The surface unit normal vector is represented by (Fig. 2.2.1) 


n 


f 


e T xe 

~i ~ 


II 


COS lb 

C 

sirn|> c 

0 


(2.3.3) 


The angular velocities of rotation are represented by (Fig. 2.2.1) 


( 1 ) _ ( 1 ) 
0 ) - 0 ) 

0 

0 

( 2 ) _ ( 2 ) 
0 ) - 0 ) 

0 “ 
0 


_ -1 _ 


1 _ 


(2.3.9) 


Vectors V^^ and 

JT 

transfer motion, 


v£^ represent the velocities of the contact point in 
while gears rotate. Thus (Fig. 2.2.1) 



i ' 


1 " 

vJV = W (1) 0,1 

-tr i 

0 

..(2) _ (2) n , 
Ytr " “ °2 Z 

0 


_ 0 _ 


0 


(2.3.10) 


Vector C and 



are 


C 


° 1°2 


0 

1 

0 


r 


( 1 ) 




= °1 I 



(2.3.11) 
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N, N N. N, 

where 0,0.. = C°=^- + -^— , 0.1=-^- 0..I = p is diametral 

12 zp 2p 1 2p 2 2p 

pitch and C° is ideal or nominal central distance. 

and are the first and the second principal curvatures of the 

gear tooth surface represented by 


(2) _ 


0 ^Ism^ 


k (2) = 0 

<H 


(2.3.12) 


W 2 

where 0^1 = is the radius of the gear pitch circle. The positive 

sign of the curvature indicates tnat the curvature center lies on tne 

positive direction of the surface unit normal. The unit vector k 2 in 

( 2 ) 

Eq. (2.3.5) coincides with the unit vector of u . Thus k 2 = k^ where 
k^ is the unit vector of the Zf-axis (Fig. 2.2.1). 

Equations (2.3.1) - (2.3.12) yield the following relation between 
tne derivative m' = (m 01 ( $. ) ) and the principal curvatures of the 

Z 1 Q <p ^ 1 J- 

crowned pinion tooth surface 


m 


21 


A 

B 


wnere 


A = 


.( 1 ) 




N^sm i/» c 


B = 


N. 


N^tanij; 


C L 


N ( 

2p + j 
Ni sln ^ c 


(1) 


2p 


- N i s mi|> c ) 


V N 2 


(2.3.13) 


The goal is to provide a parabolic type of the transmission error 
function at least in the neighborhood of the main contact point I. 
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This requirement can be satisfied with the negative value of m 2i‘ 
Using Eq. (2.3.13), we can control the value of (it is negative) 


by varying the values of k^'^. It is evident that m is negative if 


.( 1 ) 




Here : 


J2 is the curvature at point I if the 

N.sin* 

.( 1 ) 4 


N , s i n 
1 v c 

pinion would be provided with a regular involute curve; Kj X/ is the 
curvature of the modified shape of the pinion tooth in its middle 
cross-section. We have emphasized that the pinion tooth surface must 
oe deviated in the longitudinal direction (Fig. 1.3) to provide the 
point contact of the mating tooth surfaces and localize the bearing 
contact. 
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3. Design of the Crowned Pinion Tooth Surface for the 
Prescribed Transmission Error Function 

3.1 Basic Concept 

Consider that the function ^ 2 < ^ 1 ) t * iat relates the angles of 
rotation of the gears is represented by (Fig. 3.3.1). 


^1 2 
• ' 4> t ) - <l> ]_ + ^ i — a i 'J 1 


l v 1 


(3.1.1) 


N 


1 . 


Here: <p ^ — is the linear part of function 

^ 2 

<f> 2 ( (j>i) and A 2 ^ 4 > 1 ) = ^ a l^l Prescribed parabolic function of 

transmission errors. 

The gear is provided with a regular involute surface. If the 
gears are aligned, the crowned pinion tooth surface and the gear tooth 
surface are in mesh in the middle cross-section. The contact of the 
gear tooth surface can be considered as the contact of two planar 
shapes. Our goal is to determine the shape of the pinion tooth 
surface in the middle cross-section. 

We consider that the transmission errors in the region of 
- jj— < (j>^< are on average equal to zero. Thus 


It 

N 1 2 
/ jj (b i -a ]_<!>][ )d<(» 1 = 0 


(3.1.2) 


Equation (3.1.2) yields that 



(3.1.3) 
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Henceforth we will use the following coordinate systems 
(Fig. 3.1.2): s l( x l, Y l) and S 2^ X 2, Y 2^ t^ at are rigidly connected to 

the pinion and gear respectively; and Sf(xf,yf) that is the fixed 
coordinate system employed in the housing. 

The determination of the pinion tooth shape is based on the 
following considerations: 

(i) The gear tooth shape e 2 is a regular involute curve 
represented in coordinate system S 2 by the vector-function r 2 (<jb G ); 

the shape unit normal is represented by n ^ ( 4> G ) • 

(ii) Using the coordinate transformation in transition from S 2 
to Sf, we can represent the family of curves Eq in Sf as follows 

[r£ 2) ] = [M f 2 ] [r 2 ] [n< 2) ] = [ L f 2 ] [n 2 ] (3.1.4) 


(iii) The equation of meshing is represented by 


n( 2 >. V ( 21 > = 0 

~r ~r 


(3.1.5) 


( 21 ) 

where Vl ' is the sliding velocity at the contact point. 

~ t 

Equations (3.1.4) and (3.1.5) yield the following relation 


f(<f>Q/<J> 2 ) “ 0 


(3.1.6) 


that is called the equation of meshing, 

(iv) The following equation 


(rj = [M lf ][r^ 2) ] = [M lf ] [M f2 ] [r 2 ] 


(3.1.7) 
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and the equation of meshing (3.1.6) determine the shape of the pinion 
tooth in the middle cross-section. The surface of the pinion can be 
determined based on the shape in its middle cross-section. 

3.2 Derivation of Pinion Tooth Surface 
Gear Tooth Involute Shape Equations 

Tiie gear involute shape and its unit normal are represented in S 2 
as follows: 


x z = r 2 [sin<fr G -<|> G cosi/> c cos( 4> G “t|> c ) ] 


y 2 = r 2 [cos^+^cos^sint^-^)] 


(3.2.1) 


z 


2 


0 


n. 


-cos( <f> G -'J> c ) 

sinU G “* c ) 

0 


(3.2.2) 


Here: is the radius of the gear pitch circle, ^ is the pressure 

angle and <f>^ is the involute curve parameter; n£ is the normal of 
tootn shape. 

Matrix l M f 2-1 is represented by 


[M 


f 2‘ 


-cos<()2 sin^ 0 

-sincj)^ -cosrp 2 0 

0 0 1 

0 0 0 


ol 


0 

1 


(3.2.3) 
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usiny Eqs . (3.1.4), (3.2.1), (3.2.2) and (3.2.3), we represent the 

family of shapes ^ an< 3 the shape unit normal in Sf as follows 


Lr 



-sin( <t> G -<t> 2 ) + <t> G cost|) c cos( <f> G — i|> c — 4> 2 ) 
—cos ( 4> g “ 4> 2 ) ~ <f» G cosi|) sin( -<(> 2 ) 


0 

1 


(3.2.4) 


[n 



cos U G -* c - + 2 ) 

-sin( <t> G -rP c -<l > 2 ) 
0 


(3.2.5) 


Equation of Meshing 

( 21 ) 

The sliding velocity V v ' is represented by the following 
equation 


v l21) . y 2 > X (V 1 + r‘ 2 >) - e <»*r< 2 > 


(3.2.6) 


Here (Fig. 3.1.2) 
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(1) _ (1) 
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II 
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1 
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i 
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N,+N 0 

C° = r ( — 2 ) 

c 2 ' N n ; 


N 1 

“2 t0 l ( N 2 ) " 2a l*l 


The expression for is obtained by differentiating Eq. (3.1.1) 
Equation (3.2.6) yields 



r -y (2) i 

f 


” C° ' 

v (21) = 

y ( 2 ) 
X f 

, (2) (1). 

( W “U> ) + 

0 


0 


- 0 


( 2 ) 


(3.2.7) 


Equations (3.1.5), (3.2.7), (3.2.4) and (3.2.5) yield the following 

equation of meshing 

N 1 N 1 

f ( <f> 1 )= cos 'I' c (n“ + 1 “ 2a 1 * 1 ) - (-jj- + 1) cos ( <J> G — < !* 2 ) = 0 

(3.2.8) 


b- 


1 2 
where $2 = NT 'h + (b i _a i 4>i) • 


Fixiny <f> ^ , we can obtain parameter <f> G . Equation (3.2.8) provides 
two solutions for $ . Therefore it is necessary to choose the 
solution corresponding to the working part of shape E 2 . 


Pinion Tooth, Shape Equations 

Matrix [ M i f ] is represented by (Fig. 3.1.2) 


(H u i = 


COS 

sin ^ 

0 

0 


cos <|>i 

0 

0 


0 

0 

1 

0 


(3.2.9) 
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Equations (3.2.9), (3.1.7), (3.2.4) and (3.2.8) yield the following 

equations for the pinion shape 


= -C 0 sin<t> 1 + r 2 sin ) ] (3.2.10) 

N 1 2 

+ r 2 [~ 4 > 1 + b 1 - a 1 * 1 + ii> c - gUj)] cos^cos [ ( d^+g ( <\> l ) ] 

= C^cos^ - r 2 cos U 1 “'P c + g ( ^ -j. ) ] 

N l 2 

+ r 2 [|p cf> 1 + b x - a 1 4 > 1 + ^ g(<j> 1 )] cos^sin [* 1 +g ( i> x ) 1 


Here : 


2 a i<f>i 

g( 4 >1 ) = arc cos{[l- ( N y+j ] cos ^ c } 

Equations (3.2.10) represent the pinion shape in parametric form, with 
the parameter <*> 1 . For the case where the prescribed transmission 
errors are zero, a^ = 0, b^ = 0, g(<|>^) = t!> c and Eq. (3.2.10) represent 
the pinion tooth shape as a regular involute shape. 

Pinion Tooth Surface Equations 

The crowned pinion tooth surface must be deviated from a 
cylindrical surface in the longitidunal direction (Fig. 1.3). The 
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simplest way to obtain such deviation is based on representation of 
the pinion tooth surface as a surface of revolution. Such a surface 
can be generated oy tne rotation of the pinion tooth shape about a 
fixed axis. The determination of the pinion surface of revolution can 
be performed oy using the following procedure: 

(i) Assume that the pinion tooth shape has been transmitted from 
tne coordinate system to the auxiliary coordinate system S a (Fig. 
3.2.1a). Thus the tootn shape is now rigidly connected to coordinate 
system S a and is represented in S a by the following matrix equation: 


[r a ] - 


[M al ] [r l ( *l )] 


where 


(3.2.11) 


[M 


al j 


1 0 0 

0 1 U 0 

0 0 1 0 

0 U 0 1 


(3.2.12) 


(ii) On tne second stage we use a fixed coordinate and rotate 
the coordinate system S a with the tooth shape about the y^-axis. The 
tooth shape will generate a surface of revolution represented in as 
follows (Fig. 3.2.1b) 


[r b ] = 


t M ba ] [r al 


(3.2.13) 


where 
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[M 



cos 0 0 

0 1 

sine 0 
0 0 


-sine 0 
0 0 

cos e 0 

o l 


(3.2.14) 


(iii) Now, we can represent the generated pinion tooth surface 
in coordinate system S-^ using the following matrix equation (Fig. 
3.21c) : 


[ri] = [M lb ][r b ] 


(3.2.15) 


where 


[M 


lb ] 


10 0 -R. 

0 10 0 

0 0 1 0 

0 0 0 1 


(3.2.16) 


The final expression of 
represented now as follows 


the pinion tooth surface can be 


i r 1 ( 0 '*p )J 


f 1 (e, ♦p) 

£ 2 (9 'V 

£ 3 ( 0 ,^) 

1 


[Mlb] tw ba ( e) J t M al J <t> p ) 1 

(3.2.17) 


Equations (3.2.11) - (3.2.17) yield 
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f 1 ( 0 / <frp ) = X 1 {* p )CO80 “ R(l-COS9) 

f 2 (0 'V = y^d-p) 

t" 3 ( 9 / <(> p ) = x 1 ( < ()p)sin0 + Rsin9 (3.2.18) 

here x ^(<j> ), y-^(<j>p) are functions (Eq. 3.2.10) that represent the 
pinion shape in the middle cross-section. Also for convenience, we 
use <j>p as curve parameter instead of <f>^. 

3.3 Simulation of Meshing and Bearing Contact 

The authors have developed the TCA (Tooth Contact Analysis) 

program to simulate the meshing and bearing contact of misaligned 
gears (see the Appendix). The results of the TCA program confirm that 
the proposed method of crowning reduces the sensitivity of gears to 
the misalignment (as shown in numerical example). 

Example 

Given: numbers of teeth: = 20, Uo = 40; diametral pitch 

P = lu Ypj- ; pressure angle \p Q - 20°. The pinion tooth surface has been 
designed as a crowned surface and the function of transmission errors 
for the gears without misalignments has been represented as a 

paracolic function with the "Level of Kinematical Errors ( LKE ) " = 2 

arc second. Radius of rotation R^ = 350 in (see Fig. 3.2.1). The 

developed TCA program has been applied for the evaluation of trans- 
mission errors for the following misalignments: 

(i) The change of the center distance is — — = 1%. The gear axes are 

not parallel but crossed and the twist angle is 5 arc minutes. The 
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function of kinematic errors caused by the misalignments mentioned 
above is of a parabolic type and the maximum value of transmission 

i 

j. errors is 1.2 arc seconds. 

(ii) The year axes are not parallel but intersected and form the 

I 

angle a = 5 arc minutes. The function of kinematic errors is of a 
parabolic type and the maximum value of transmission errors is 2.0 arc 
seconds. 


I 


i 

i 

i 

j 

i 

i 

i 


i 

I 

i 


l 
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4. Deviation of Pinion Tooth Surface by Tool with Cone or 

Revolute Surface 

4.1 Introduction 

A new method for generatation of the pinion of spur gearing is 
developed. The involute spur gears with the crowned pinion tooth 
surface are less sensitive to the gear misalignment, the bearing 
contact has a favorable location and the kinematical errors are very 
small and witn favorable shape. The crowning may be done by grinding 
or shaving using a tool with conic surface, a surface of revolution or 
a planar surface. The new method covers: (i) the theory of the 

proposed method of crowning, (ii) determination of the pinion surface 
generated oy the tool, (iii) conditions of non-undercutting, (iv) 
principal curvatures and directions of the pinion tooth surface, (v) 
dimensions and orientation of the instantaneous contact ellipse, 
(vi) the bearing contact and (vii) the kinematical errors due to 
gear misalignment. Computer programs are developed in all stages 
including the simulation of meshing and bearing contact TCA (Tooth 
Contact Analysis) program. 


4.2 Principle of Generation and Coordinate Systems 

Consider two rigidly connect generating surfaces Zq and E^. The 

generating surface is a plane and generates the gear tooth surface 

E 2 that is a regular involute surface. Surface E^ is a surface of 

revolution. Initially we consider that E is a cone surface and 

P 

Ep and E^ contact each other along a straight line that is the 
generatrix of the cone. 

Figure 4.2.0 shows the generating surfaces E p and z . Figure 1.4 

o p 
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snows a sKetch of the tool that could he used for pinion generation. 

I figure 4.2.1 illustrates the process of generation. While the rigidly 
connected generating surfaces perform a translational motion, the 
pinion and gear rotate about their axes 0^ and 0 2 respectively. The 
parameter of motion of the cutters, s, and the angles of rotation of 
the pinion and the gear, <f>^ and <f> G , are related as follows 

l 

S = r l*p = r 2*G (4.2.1) 

| where r^ and r 2 are the gear centrodes radii; the cutter centrode is 
the straight line that is tangent to the gear centrodes. Point I is 
the instantaneous center of rotation. Coordinate systems and S 2 
are rigidly connected to the pinion and the gear respectively. The 
riyicily connected tool surfaces (see Fig. 4.2.0) are represented in 
coordinate system S c ; Sf is the fixed coordinate system. 

The generating surface E G (a plane) is covered with the set of 
contact lines L G2 — the instantaneous lines of tangency of surfaces 

| Eq and E 2 (Fig. 4.2.2a). The location of these lines depends on the 

j value of parameter Line L Gp is the line of tangency of generating 

! 

surfaces E Q and E p . Points , M 2 ,..., M n represent the set of points 
of intersection of L Gp and L G2 . The coordinate transformation from S c 
to S 2 is represented by the equation 

[r 2 ] = [M 2c ] [r^ G) ] (4.2.2) 

n(G).v(G2) = £ 2 (u G ,d G'*G ) = 0 (4.2.3) 
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Here: ^ U G , 0 G^ the vector function that represents plane 

e G ' ^ M 2 c j is the matrix that describes the coordinate transformation 
from S Q to 82 * Equation (4.2.3) is called the equation of meshing of 
the gear cutter and gear 2 [Litvin, 1987]. n v 'is the unit normal of 
generating surface and is the relative velocity 

/v /v /v ■* 

for the process of generation. Contact point Mj_ of surfaces and 
1 2 may be represented in coordinate system S c as a function of 
parameter <j> G as follows 


£c ^ £c ^ U G ,0 G^ f 2^ u G ,9 G'*G^ " 0 

F 2 (Ug/0g) = 0 


(4.2.4) 


Equations 


r (G) 

~C 

= ^ G)(u g' 9 g>' 

F 2 (U G ' 0 

G> = 0 



(4.2.5) 

represent 

line Lqp in 

coordinate 

system 

S C‘ 

Equations 

(4.2.4) 

represent 

the location 

of points M 

1 , M 2 ,. 

. . , M n 

on surface 

Z G as a 

function 

of parameter 

4>G • Matrix Eq. 

(4.2. 

2) and Eq. 

(4.2.4) 

represent 

the set of contact points 

of gear 

tooth 

surfaces 


Z 2 and Z^ 

in coordinate 

system S 2 « 

These 

points 

are the cen 

ters of 


the instantaneous contact ellipses that form the bearing contact of 
gear 2 and pinion 1 . 

Using the matrix equation 

[r f ] = [M fc ] [r c b)] (4.2.6) 
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and Eq. (4.2.4) we may represent in coordinate system Sf the line of 

action the set of contact points of E and E^ in the fixed 

coordinate system. 

Similarly, considering the generating surface E^, we may 
represent in coordinate system S c the lines of contact, Lp^, of 
surfaces E p and E as follows 



Ec P)(u p' 9 p ) 


n^ P ^.v^ P1 ^ = f 1 (u p ,Q p , 


,) = 0 


(4.2.7) 


( P ) 

Here n v ' is the unit normal of generating surface 

E„; V (P1) = V^ p ) - is the relative velocity for the process of 

generation . 

Tne pinion surface is represented by the equations 


up = 


‘"lc* 


-‘ P> > - 


[M lf ]|Mf c ][r‘ p >] 


f^ (^p, 0p , 


,) = 0 


(4.2.8) 


Equations 

£c P) = Sc^^P'^' p ( u p' 9 p ) “ 0 (4.2.9) 

represent, in terms of u p , and e p , the line of contact of surfaces 
Ep and E g » in coordinate system S c . Contact point of surfaces 
E ^ and E 2 raay be represented in coordinte system S c as a function of 
parameter <j> p as follows 

P P) = P P) ( Up , 9 p ) fpup^p^p) = 0 (4.2.10) 


F 1 (Up,0 p ; 


0 
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Using the matrix equation 


[r f j = [M fc ] [rJ, P) ] (4.2.11) 

and Eq. (4.2.10) we may represent the line of action of gears 1 and 2 
as a function of <j>p in coordinate system Sj. Equations (4.2.4) and 
(4.2.6), and (4.2.10) and (4.2.11), determine the same line of action 
but in different parameter, <f>g or <f>p, respectively. 

4.3 Tool Surface 

The pinion generating surface is a cone and may be represented in 
an auxiliary coordinate system S a as follows (Fig. 4.3.1) 

x = u n sinacose n y = d-u n cosa (4.3.1) 

a P P J a P 

z a = UpSinasinQp 0<u p < ^ 0 <e p <2n 


The surface normal is represented by 


3 r, 

~a = Tel 


3r 

ts; = u p sina 


cosaCOSQp 

sina 

cosasin9p 


(4.3.2) 


The unit surface normal is (provided Upsina * 0) 


Ln a ] 


cosa cos9p 
sina 

cosa sin0p 


(4.3.3) 
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Figure 4.3.2 illustrates the installment of the tool cone in 
coordinate system S c . Axes of coordinate system are parallel to 
axes of system S a and axes Xj.,, y^ lie in plane z c = 0. The coordinate 
transformation from S a to S c is represented by the following matrix 
equation 




t»ob> lM tJ lr a> = ‘"ca^a 1 


(4.3.4) 


Here 


IM cb> - 


cos ( a-* ) 

V 

■sin ( a-i^ c ) 
0 
0 


sin(a-']; c ) 
cos ( a-6 ) 
0 
0 


0 

0 

1 

0 


0 

0 

0 

1 


(4.3.5) 


[M ba ] " 


1 

0 

0 

0 


0 

1 

0 

0 


0 

0 

1 

0 


-UpSina 

* 

-(d-UpCOSa) 

0 

1 


(4.3.6) 


Combining Eq. (4.3.5) and Eq. (4.3.6) results in: 


[M J = 
ca 


cos (a-tp c ) 
■sin (a-<fr c ) 
0 
0 


sin( a-<P c ) 
cos ( ct-i/i ) 
0 
0 


0 

0 

1 

0 


-UpSin^-dsinf cx— ) 
* 

u pCos i^-dcos ( a- ij> c ) 
0 


(4.3.7) 
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d_ 

cos a 


* 

Here u p 



b 

COS ’i 
r C 


Using Eqs. (4.3.4), (4.3.7) and (4.3.1), we represent 

the generating surface E p in coordinate system S c by equations as 
follows 


[r 



* 

u p [sinacos9 p cos(a-i|) c ) - cosasinfa-^)] - u p simji c 

* 

-u p [sinacos0 p sin(a-i|» c ) + cosacosta-^)] + u p cosi|> c 
u p sinasin9 p 

1 


(4.3.8) 


It is evident that the point of the cone surface with the surface 

* 

parameter u p = u p and e p = 0 coincides with the origin 0 C of 
coordinate system S c . The unit normal of cone surface is represented 
in coordinate system S c by 


[n 



[L j [n 

ca a 


(4.3.9) 


Here 


[L 


ca‘ 


cos ( a-ip c ) 
-sin ( a-ij» c ) 
0 


s in ( a-i|/ c ) 0 

cos(a-i); ) 0 

0 1 


(4.3.10) 


Equations (4.3.9), (4.3.10) and (4.3.3) yield 
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[n 



cosacos 0 p cos ) + sin a sin (a~\l> c ) 

-co3acos0 p sin( a -ii) c ) + sin a cos (a“i!> c ) 
cosasin 9 p 


(4.3.11) 


4. 4 Pinion Tooth Surface 

The equation of meshing of the generating surface and the 
pinion tooth surface is represented by 


n (P). V (PD = n< p > 




(4.4.1) 


We may also use 
contact normal must i 
(Litvin, 1987). Thus 


X -x (P) 
c c 


Y -y (P) 
c 1 c 


n 


(P) 


cx 


n 


(P) 

'ey 


the equation 
ntersect the i 
we obtain (Fig 



-z (P) 
c c 


n 


(P) 


cz 


based on the condition that the 
nstantaneous axis of rotation I-I 
. 4.2.1) 


(4.4.1a) 


Here: (X x , Y c , Z c ) are the coordinates of a point that lies on axis 

I-I; X' , Y' and Z' are the coordinates of the cone surface; n ov , 

c C C 

n C y and n cz are the projections of the surface unit normal. Here 
(Fig. 4.2.1) 


X c = 5 = r^p Y c = 0 

Equations (4.4.1a) (4.4.2), (4.3.8) and (4.3.11) yield 


(4.4.2) 
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f l ( u p» 9p» <f>p) = r 1 'f’p i “cos a cos e^sin ( a-i() c ) + sinacos ( ot— if> c ) ] 


* 2 2 
■U D COS0 D + u u (cos9 d cos a + sin a) = 0 


(4.4.3) 


Here (Fig. 4.3.2) 


* 



AO 

c 


_d b_ 

cosa cost|) c 


(4.4.4) 


The line of contact LpQ of generating surfaces E p and Eq 
(F ig. 4.2.0) is generatrix of the cone surface determined with 


9 p = 0 (4.4.5) 

Equations (4.3.8), (4.4.3) and (4.4.5) represent the location of the 

instantaneous contact point of surfaces E^ and E 2 on surface E p as 
follows 

X c = r i^p sin * c y c = ~ r^psm^cos^ = 0 

(4.4.6) 

The coordinate transformation in transition from S c to Sf is 
represented by the matrix equation 

[r f j = [M fc ] [r£ P) ] (4.4.7) 

where (Fig. 4. 2.1) 
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(4.4.3) 


t M fc ] 


10 0 -r^p 

0 10 r.^ 

0 0 1 0 

0 0 0 1 


Equations (4.4.b) to (4.4.8) yield 


x 


f 


(|>pCOS 



y f = r 1 " r i<l , pSi n ’J’ c cos<j, c z f = 0 (4.4.9) 


Equations (4.4.9) represent the line of action as a tangent to the 
base cylinder of radius r^ that lies in plane Zf = 0 and passes 

through point I (Fig. 4.4.1). 

The pinion tooth surface may be represented in coordinate system 
S-^ as the set of contact lines of surfaces Ep and E-^. Thus we obtain 


irp - W lc ][r‘ P) ] - [M 1 £ nMf c ][r‘ p >] 

f i ( u p/ 0 p » «J>p) = 0 


(4.4.10) 


Here (Fiy. 4.2.1) 


[Mif] 


cos^p sin<{>p 0 0 

-sinifip cos cf>p 0 0 

0 0 10 

0 0 0 1 


(4.4.11) 


Equations (4.4.10), (4.4.3) and (4.3.8) yield 
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X 1 = u p [ cos Qpsinacos ( a-ij> c + <i>p) - cosasin ( a-^ c +4 p ) ] 

* 

-u p sin( <p c -<j> p ) + r 1 (-<| )p cos<|> p +sin<f, p ) 

y-j_ = - u p [cose p sinasin( a-^ c +<() p ) + cosacos ( a-i|> c + <{> p ) ] 

* 

+ Up cosU c -$ p ) + r L ( <J>. p sin<f> p + cos<f> p ) 

Z 1 = u p s: '- n( ) , p s i na 

f i (u p/ 9 p , <(»p ) = r i 4> p [-cosacos 0 p s in ( a- ip c ) + sinacos ( a- 1 |> )] 

* 2 2 ( 4 4 1 o \ 

-UpCOS0 p + u p ( cos 0 p cos a + sin a) = 0 v 1 *****-*-^ 

Equations (4.4.12) are equivalent to the representation of the pinion 
surface by the vector function r ^ (u p , 0 p , <f> p ) and the equation 
of meshing f ^ ( u p , 9 p , <j> p ) = 0. Taking it into account that 
f i ( u p ( 0p» <J>p) = U is linear with respect to the parameter u p , it is 
easy to eliminate u p and represent the pinion surface in two 
parametric form, with the parameters, 0 p and <(i p . 

The intersection of the surface by the plane = 0 represents a 

regular involute curve. Equations of this curve may be derived from 
the surface Eq. (4.4.12), taking it into account that 

★ 

0 p = 0 and u p -u p = r ^ <j> p s i n (4.4.13) 

Thus, we have 

X 1 = r l t-4’p cos ^ c cos( + sin< frpl (4.4.14) 

y l = r i U p eos* c sin( + cos^ p ] 
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Cutting the pinion tooth surface by a plane z-^ = const, we obtain the 
cross-section of the pinion tooth surface that is deviated from a 
regular involute curve (Fig. 4.4.2). 

It will be shown later that the kinematic errors caused by the 
gear misalignment can be of two types that are shown in Fig. 1.2a and 
1.2b. Kinematic errors shown in Fig. 1.2a correspond to the case when 
tne transformation of rotation is interrupted and the gear tooth 
surfaces are out of contact. To avoid the appearance of such 
kinematic errors, a surface of revolution that slightly deviates from 
the cone surface is actually used. The determination of curvatures of 
the surface of revolution is a subject of optimization. Since a 
surface of revolution is used instead of a cone surface, the crowned 
pinion surface deviates from a regular involute surface in all cross- 
sections including z^ = 0, the middle cross-section of the tooth 
surface . 

4.5. Conditions of Pinion Without Undercutting 
General Approacft 

The proolem of undercutting of the pinion tooth surface by 
crowning is related with the appearance on the pinion tooth surface of 
singular points. It is known from differential geometry that the 
surface point is singular if the surface normal is equal to zero at 
such a point. Litvin, proposed a method (Litvin, 1987) to determine a 
line on the tool surface which will generate singular points on the 
surface that is generated by the tool. This line designated by L 
(Fig. 4.5.1) must be out of the working part of the tool surface to 
avoid undercutting of the pinion by crowning. 
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The limiting L of the tool surface is determined by the following 
equations 




9 P ) 


(4.5.1) 




(4.5.2) 


F(u p/ e p ,$ p ) = 0 


(4.5.3) 


Vector Eq. (4.5.1) represents the tool surface, Eq. (4.5.2) is the 


equation of meshing (see Eq. (4.4 
requirement that (Litvin, 1937) 


au p 

i x i P) 

a 9 p 

_v (cl) 

cx 


Jp) 

9y c 

^ p > 

_ v ( cl > 


9U P 

30p 

V cy 


3f l 

9f l 

3£ 1 

d^>p 

a Up 

a 0 p 

a ^p 

dt 

9y c 

*4 P) 

-v (cl) 

cy 


a u p 

a 0 p 


■«‘ P> 

32 c P) 

-v (ci) 

cz 


9U P 

3 9 p 


9f l 

9f l 

_ !£i 

d if) p 

9U P 

a 0 p 

a <Pp 

dt" 


3)) and Eq. (4.5.3) follows from the 


»x< P) 

-< p) 

_ v ( cl ) 

CX 


aup 

a 0p 


»‘ P) 

92 < P) 

-v (cl) 

CZ 


9 Up 

a 0p 


9f i 

9f i 

9f l 

u <j> p 

au'p 

a 0 p 

a<f> P * 

dt 


(4.5.4) 
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Derivation of Equation 


( Cl ) 

Relative velocity V c is to be represented in coordinate system 


b c as : 


,(cl) 


= V (c) - v (1) 


(4.5.5) 


wnere is the velocity of the cutter and ^ is the velocity of 

~c ~c 

the pinion at generating point. Using the sketch in Fig. 4.2.1, we 
get 


V" 

~c 



- ds 
dt 


i 

a* 

3 

rH 

1 

1 

= 

0 

= 

0 


1 

o 

i 


1 

o 

1 


(4.5.6) 


r (D 


= K P x £ c + °c °1 X ~P 


(4.5.7) 


Deriving Eq. (4.5.7) we substituted the sliding vector nip that passes 
through U-^ by the equal vector tnat passes through 0 C and the vector- 
inome nt U c 0 ^ x w p . 

Here : 



■ o ■ 


r i*p' 

— 

0 

II 

jd 

|o 

3 

-r n 



C 1 

1 


1 


0 


Equations from (4.5.6) to (4.5.8) yield 
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(4.5.9) 

where ^ and y£ P ^ are represented by equations (4.3.3). Equations 
(4.5.9), (4.3.3), (4.4.3) and (4.5.4) yield 

2 

F l(u p /0p) = m 1 m 2 q - (m 2 +m 1 m 4 )qp - m 3 m 5 q+m 4 p = 0 (4.5.10) 


V 


(Cl) _ 


-x 


if) 

c 

(P) 


+r l*p 


Here 


Q = 


UpCOSa 


P = 


U n COS 

P i_c 
d/cosa 


m 


4 = sinasin ( a-^ c ) cos e p + cosacos ( a~'P c ) 


= sinacos ( a-ip ) - cosasin ( a-i|i ) cos J 0 p 


m 3 = [sinacos ( a-i|) c ) - cosasin ( a-i|> c ) cos 9 p ] 


• 2 

m 4 = sinacosasin 0 p 


r ^cosa 


m 5 <f 


8F- 


Using Eq. (4.5.10) and taking it into account that , . * 0 , we 

d 0p 

can represent u p as function (0p). Undercutting will be avoided if 
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(4.5.11) 


vv 



cos a 


The analysis of Eq. (4.5.10) yields that inequality (4.5.11) is 
satisfied if the conditions of non-undercutting of the pinion 
generated by a regular rack cutter are satisfied. 

4.6. Principal Directions and Curvatures of Tooth Surfaces 
Introduction to Pinion Principal Directions and Curvatures 

The direct determination of pinion principal directions and 
curvatures requires complicated derivations. A simplified approach 
for the solution to this problem has been proposed by Litvin, in 
19 69 . The main idea is as follows: 

Consider that the tool surface E p generates surface E^ as the 
envelope of the family of surfaces E p . Let Lp^ be the instantaneous 
line of contact of E p and E^ and M is a point of this line. We 
consider as given the principal curvatures and directions of E p at any 
point and also the parameters of motion of surfaces E p and E being in 
mesh. Then the principal curvatures and directions of the generated 
surface E can be expressed in terms of principal curvatures and 
directions of the generating surface E p by using the equations that 
have been derived by Litvin, 1969. 


Principal Curvatures and Directions of Ep 

The tool surface is a cone surface and its principal directions 
coincide with the direction of the cone generatrix and the direction 
that is perpendicular to the cone generatrix. The Rodrigues' formula 
(Eq. 4.6.1) can be used to find the principal curvatures and 
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directions . 


K 


i,n*~r 



(4.6.1) 


Here: < T tt are the principal curvatures of surface. V is the velo- 
1 , 1 1 ~r 

city of a point that moves over a surface and h is the derivative of 
the surface unit normal n, when n changes its direction due to the 
motion over the surface. 

Equation (4.6.1) yields the following expressions for the 
principal curvatures and directions (see Eq. (4.3.1) for the tool cone 
surface S p represented in coordinate system S a ) 


-sine p 
0 

cose p 

K (P> 1 

I u p tana 

cos 9 p sina 
-cosa 
sine p sina 



(P) _ . 9 _£p 

~ii 3 u p ’ au p 


(P) _ 

'I 3 6 t-j 


3r . 


(4.6.2) 


(4.6.3) 


(4.6.4) 


(4.6.5) 


The negative sign for indicates that the curvature center is 

located on the negative direction of the surface normal. 

The principal curvatures are invariants with respect to the used 
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coordinate system. The unit vectors of the principal directions, 
e j. ' amd , may be represented in coordinate system Sf by using the 

following matrix equation 




(4.6.6) 


( p) / p \ 

Here: le' '] and [ei. ' ] ^ are vectors represented in S and S £ , 

respectively; matrix [L ca ] is represented by Eq. (4.3.10); matrix 
LL,£ C J is the 3x3 unitary matrix (see matrix Eq. (4.4.8)). 

Equation (4.6.6) yields 


Le 



f 


-sine p cos ( a-\l> c ) 
sin0 p sin(a-^ c ) 

COS0 p 


[e (P) ] 
ie II J f 


cos0 D sinacos(a-i|» ) - cosasin( a-rt ) 
-cos0 sinasin(a-if) ) - cos«cos(«-* ) 

Jl 

sin0 p sina 


(4.6.7) 


(4.6.8) 


Principal Curvatures and Directions of Pinion Tooth Surface 

The determination of principal curvatures and directions for the 
pinion tooth surface is based on the following equations (see Litvin, 
iyu7) 


tan 2 ct 


(Pi) _ 


2b 13 b 23 


.2 ,2 , (P) (P) v. 

b 23 b 13 ( K I K II b 33 


(4.6.9) 
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(4.6.10) 


( 1 ) 

II 


,2 2 (P)_ (P) v h 

(1) _ b 23 b 13 ( K I K II b 33 


b 3 3 c °s 2 a 


(PI) 


(1) ^ (1) _ (P) ^ (P) . b 13 +b 23 

< rT + k t " < T +k tt + K 

XX I I II D33 


(4.6.11) 


( P ) ( P ) 

Here: and Kjj are the principal curvatures of the tool surface 

( P ) ( P ) 

E n ; e' and e' T ' are the unit vectors of the principal directions on 
(Fig. 4.6.1); and e| b and eii^ are the principal 

ir 1 II ~ 1 1. 

( Pi ) . 

curvatures and directions on the pinion tooth surface E ^ ; and a is 

measured counterclockwise from e|. P ^ to (Fig. 4.6.1). The 

coefficients b 13 , b 23 and b 33 have been derived from equations 

represented in [Litvin 1937 ] but modified for the case when a rack 
cutter generates a gear. The expressions for b^ 3 , b 23 and b 33 are as 
follows 


b 


13 


p ( p ) 

Sirs] 



0 


b 


23 



0 


“ K 


(P) 

II 


V 


(PI) 


r (Pl) 

II 


(4.6.12) 


b = [no, V< P) ] + [no, V (P1) ] 

J J ~ p ~ L L p 


. (P)(V (P1 )) 2 - K (P) (V (P1) ) 2 - -a 
K I v I ' K II K II ' m 


IP 


"ip bE (P) i f ' 


(4.6.13) 


where m^ p = — = 

s s 

The vectors that are used in Eqs. (4.6.12) and (4.6.13) are 
represented in coordinate system Sf (Fig. 4.2.1); i^, and k^ are 
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The expressions of 


tne unit vectors of coordinate axes of system Sf . 
tnese vectors are as follows 


w-cokj; (4.6.14) 

is the angular velocity of the pinion being in mesh with the rack 
cutter 


V 


(P) 

tr 


-a,r lif 


(4.6.15) 


is the transfer velocity of a point of the rack cutter that performs 
translational motion; r^ is the radius of the pinion centrode 



i 

j 

k 


* 

1 

K 

l-h 

1 

v (1) _ VK (P) _ 

~tr “ x £ 

0 

0 

(0 

- w 

x f 


X f 

4-1 

>i 

Z f 


0 


is the transfer velocity of the pinion. The term, "transfer velocity" 
means that the velocity of a point that is rigidly connected to the 
tooth surface is considered. Such a point performs the motion with 
the tooth surface 


„(P1) _ „(P) v (l) _ 

V = V - V -to 



(4.6.17) 


is tne "sliding" velocity — the velocity of a point of the rack cutter 
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can be 


( Pi ) 

with respect to tne same point of the pinion. Vector V 
represented in terms of components x Q and y c using matrix Eq. (4.4.8) 
that describes the coordinate transformation from S Q to Sf . Then we 
obtain 


x f = 


- r 


1<p p 


= 


Y + r , 
1 c 1 


(4.6.18) 


Equations (4.6.17) and (4.6.18) yield 


V 


(PI) _ 

- u 


^c 

_x c +r l*P 


0 


(4.6.19) 


Tne principal curvatures of the cone that is rigidly connected to the 
rack cutter are represented as 


J P) 1 

I u tanoi 


(P 

I 


0 


(4.6.20) 


The unit vectors of principal directions are designated by 

ej. P ^ and e^^. Components Vj P ^ and Vj P ^ are represented as follows 


V (PD = v < Pl >.e (P) 


V 


(PI) 

II 


V ( P1 ) 



(4.6.21) 


After suostitution we obtain 


'13 


usine, 


sina 


u, 


U P 2 r l^P 
[1 - — cos a + — — cosasin( a-A ) ] 
Up Up c 


(4.6.22) 
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(4.6.23) 


b 23 = " “ cos9 p 


*^33 = w ^i + ^2 + ^3 + ^4^ 


where 


• 2 


k^ = UpSin QpCosasina - u p sinacosa ( 1-cos e p ) 


k 2 = r i < t , pt cos0 p cos ° icos ( a-i|) c ) + sinasin( a-t|» c ) ] 


k 3 = r l E cos 0p cos ot s in ( a-ii> c ) - sinacos ( a-i(j c ) ] 


sin 2 0 p * 

K 4 = IT^T^ l(u p -u p )o°s« - r^slnfc-^)] 


For tne contact point located in tne middle section of the 
have 9 p = 0 and 


b 13 " 0 


'23 


- 0 ) 


b 33 = “ lr l*P cos + c ” r l sin * c l 


Then we obtain the following expressions for the pinion 
curvatures and directions 


a (Pl) = 0 


K d) = K (P) = 1 

I I utan^ 


(4.6.24) 

(4.6.25) 

(4.6.26) 

(4.6.27) 

(4.6.28) 
pinion we 

(4.6.29) 
principal 
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„(1> = .(PI 1 = , -1 

II II r i sl n* c -^ 1 +pC° 3 » c . r 1 sint c -i 1 <tp°°s 






U) = e (P) 

II ~I 


(4.6.30) 


The principle curvature is the same as the principal curvature 

of a regular involute surface. However differs from zero 

because the pinion is generated by a cone but not a plane as in the 
case of a regular involute surface. 


Gear Principal Curvatures and Directions 

The determination of the gear principal curvatures is based on 
the relation between the principal curvatures for spur involute 
gears. Figure 4.6.2 shows that two mating involute shapes are in 
tangency at point M. Let us designate the principal curvature for the 
pinion by and where (see Eq. 4.6.29) 


K (U 1 K d) L__ 

I utan^ c II r 1 sini|) c -r (fipcos i|i c 


(4.6.31) 


The negative sign of and means that the radii of curvatures 

are opposite to the direction of the surface normal. The principal 

( 2 ) 

gear curvature = 0. The radii of curvature of spur involute gears 

are related as follows (Fig. 4.6.2) where <1^ = - — — 

a II r.sinti-r, ^cos ii) 

1 r c IP c 


( 1 ) 

'II 


( 2 ) 

'II 


KL = rb l+Cb2 = r l +r 2 
tan^ sin^ 


N l + N 2 
2Psin ip 


(4.6.32) 
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Then we obtain 


( 2 ) 


r l +r l 

sin*. 


K 


( 1 ) 

II 


N l +N 2 

2psin^ 


( 1 ) 

II 


(4.6.33) 


The angle a 


( 21 ) 


between the principal curvatures is equal to zero. 


4.7. Contact Ellipse and Bearing Contact 

The dimensions and orientation of the instantaneous contact 
ellipse can be determined on the basis of equations that have been 
developed by Litvin (Litvin et al. 1982, and Litvin 1987). The input 
data for the computation is: ^ (i=l,2) 


and e where 


(i) 


and 


(i) 


(12) , 


jj are the principal curvatures, a ' is the angle formed 

by the unit vectors of principal directions, e^^ and e| 2 ^ and e is 

the elastic deformation of the contacting surfaces at contact point. 

The bearing contact is formed by the set of instantaneous contact 

ellipses that move over the gear tooth surfaces in the process of 

meshing. In our case, when the crowned pinion is meshing with a gear 

without misalignment, ej. 1 ) = ej. 2 ^ an< 3 cr^ 12 ^ = 0 (Fig. 4.7.1). The 

axes of the contact ellipse are directed along the n~axis and 

g-axis , respectively. Generally the orientation of the contact 

ellipse is determined with the angle a^. In the case where 
( 12 ) 

a = 0, a is equal to zero. The equations to be used for the 

computation are as follows 


= 4 u ( 1 , -k (2 > - 


9 1-9 2 1 


(4.7.1) 
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B = I [k 1 1)-k( / ) + l g l" g 2 l ] 


(4.7.2) 


where 


: (1) = <i 1} + 


II 


k < 2 > - < i 2) + 4 2) 

e I II 


(4.7.3) 


a,- h 1 ’ - kP> 


'1 I 


'll 


- J 2 ) _ ( 2 ) 

K II 


g 2 K I 


(4.7.4) 


a* - e V2 
a " A 


b B 


1/2 


(4.7.5) 


The equations for ^ and (i = l,2,) have been represented in 

Section 4.6. Figure 4.7.2 illustrates the bearing contact for spur 
gears with the crowned pinion tooth surface. 


4.6. Simulation of Meshing and Determination of Kinematical Errors 
Coordinate Systems 

The simulation of meshing is a part of the computer aided TCA 
program developed by the authors. The simulation of meshing is based 
on equations that provide the continuous tangency of contacting 
surfaces . 

Henceforth we will use the following coordinate systems: (i) Sf 
that is rigidly connected to the frame and (ii) an auxiliary 
coordinate system S^ that is also rigidly connected to the frame, and 
(iii) systems S-^ and S 2 that are rigidly connected to the pinion and 
the gear, respectively. The origins and coordinates axes and z-^ of 
coordinate systems Sf and coincide with each other (Fig. 4.8.1a). 
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' Axis Zf is the axis of rotation of the pinion. Also, the origins and 
coordinate axes z^ and Z 2 coincide with each other and z^ is the axis 
of year rotation (Fig. 4.3.1b). The errors of assembly of the gears 

i 

i 

I are simulated with the orientation and location of coordinate system 
S n with respect to Sf . Figure 4.8.2 shows the orientation 
| (a) The gear axes are crossed (Fig. 4.8.2a) 



-cosAy 0 

sinAy 

0 

lM £n J * 

U -1 

0 

C 


sinAy 0 

cos Ay 

0 


0 0 

0 

1 __ 

(b) The gear 

axes are intersected 

(Fig. 

i 

O 

i — 1 

1 

0 

0 

i = 

0 -COSA6 

-sin A6 

C 

i 

0 -sinA6 

cos A6 

0 

! 

O 

o 

0 

1 


(4.8.1) 


(4.8.2) 


Pinion Tooth Surface Equations 

T'ne pinion tooth surface has been represented in coordinate 
system by Eq. (4. 4. 12). The coordinate transformation form S-^ to Sf 
is represented as follows 

[r< 1} ] = [M fl ] [r x \ (4.8.3) 


Here 
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0 


0 


[M 



cos (f> sin<j>^ 

-sin^ cos^ 0 

0 0 1 

0 0 0 


0 

0 

1 


(4.0.4) 


where <f>-^ is the angle of rotation of the pinion being in mesh with the 
gear . 

Equation (4.4.12) represent the pinion tooth surface in a three- 
parametric form in terms of parameters u p , 9 p and <f>p where <f>p is the 
angle of rotation of the pinion being in mesh with the coned cutter. 
These parameters are related by the equation of meshing and only two 
of them are independent. It is easy to eliminate parameter u p and 
represent the pinion tooth surface in in the two-parametric form. 
Thus : 


= £iU, 


» 6p) 


(4.8.5) 


where <f>p and 0 p serve as the surface parameters. The pinion tooth 
surface unit normal can be represented in coordinate system Sf as 


In 



[ L f 1 J [n 1 J 


(4.8.6) 


The 3x3 matrix [ L,^ ^ J can be developed from [M^^J by eliminating the 
last row and the last column. The surface unit normal, n^ , is 
represented in by the following matrix equation 


tn x J = 




[L iP 1 L f C i In c P>1 


(4.8.7) 
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( p \ 

Here: is represented by Eqs . (4.3.11), [Lf c ] i- s the 3x3 unitary 

matrix and [L^f] is the 3x3 submatrix of matrix [M^f] given by Eq. 


Tne final expression for n^^ is as follows 
in* 1 ), = [L fl HL 1£ ][L fc ][n< P »] 


(4.3.8) 


Here: [Lf c ) tne 3x3 unitary matrix [I], but [Lf-^] [L^f] * [I] 

because the elements of these matrices are expressed in different 
terms - <f>p and <j>^ 

Gear Tooth Surface Equations 

The year tooth surface is a regular involute tooth surface and 
may be represented in coordinate system S 2 by equations that are 
similar to Eq. (4.4.14) 


x 2 = r 2 [-<t> G cos\p c cos ( <J>q-^ c ) + sin<j> G ] 
y 2 = r 2 [<j> G cos^ c sin( <t> G “^ c ) + cos<j> G ] 

z 2 = X (4.8.9) 


where <f> G is the angle of rotation of the gear being in mesh with the 
rack cutter. The gear tooth surface unit normal is represented by the 
equations 


n . 


W.. 



3 r ,, 3 r .. 

~2 _ X 3 £ 


(4.8.10) 


Equations (4.3.9) and (4.8.10) yield 
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(4.3.11) 


L n 2 i 


— cos ( <t> G ~t c ) 
sin ( * G -$ C ) 
0 


To represent the gear tooth surface and its unit normal in 
coordinate system we use the following matrix equations (Fig. 

4.3.2) 


|r‘ 2) i = [M h2 ] [r 2 ] 


[n< 2) ] = [L h2 ][n 2 ] 


(4.8.12) 


tM h2 ] = 


cos $2 
s i n <f> 2 


0 

0 


-sin<j> 2 

cos<ji 2 

0 

0 


0 0 

0 0 

1 c 

0 1 


(4.8.13) 


where <j> , is the angle of rotation of the gear being in mesh with the 
pinion. Equations (4.8.9), (4.8.11), (4.8.12) and (4.3.13) yield 


( 2 ) _ 


h 


= - r 9 {<(> r cos<|) cos [( <J>~-<J> r ) + )] + sin (<j>,.-<j> r ) } 


( 2 ) 

^h = 


= - r 2 { (figCos i(i sin[( <! ) 2 - () ) G ) + “ cosU 2 -0 G )} 


z (2) = A 
z h A 


(4.3.14) 


< 2> ‘ - 


-cos [ ( <t> 2 ~<P G ) + ^ c ] 

-sin [ ( <f> 2 “<f> G ) + t c ] 

U 


( (4.3.15) 
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X = 0 determine 


( 2 ) 

These equations with parameters ^ 2 ~ = ^ ' Z h = 

point M with coordinates (0, r 2 , 0) and the surface unit normal at M 
with components (-cosi/; , - siniK, 0) (Fig. 4.8.3). 

To determine the gear tooth surface equations and the surface 
unit normal in coordinate system Sg we use the following matrix 

I equations 

l 

i 

| ir^ 2) ) = LM fh j [r^ 2) ] [n< 2) ] = [L fh Hn£ 2) ] (4.8.16) 

I 

j 

Matrix [M^] (see Eqs. (4.8.1) and (4.8.2)) describes the coordinate 

I 

transformation from S h to S^.- for the case where the gear is misaligned 
with respect to the pinion. Figure 4.8.4 illustrates the case where 

' the gears are not misaligned, the shortest distance of gear axes 

C = r^ + r 2 and the tooth surfaces of the pinion and the gear contact 
each other at point M. The misalignment of the gear will not cause 
the surface contact at the tooth edge since the pinion tooth surface 
is crowned and deviates from a regular involute surface. 


Equation of Tangency 

The contact of gear tooth surfaces is simulated in the developed 
TCA program by the following equations 


.( 1 ) 


£f ( <Pp' $1 ) 



a , <p 2 ) 


(4.8.17) 


if 1 ’ 


ni ' ( <f> 


p / 9p f <J>2 ) _ 


( 2 ) 
— Y~\ 1 / 


Sf Uq'^ 


(4.8.18) 


At the contact point vector Eq. (4.8.17) provides the equality of 
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position vectors and vector Eq. (4.8.18) provides the equality of 
surface unit normals. Vector Eq. (4.8.17) provides three independent 
scalar equations that relate <f»p, 9 p , <f> 3 , <j> G r X and j> 2 * However, vector 
Eq. (4.8.18) provides only two independent scalar equations since 
n^^ and are unit vectors and | n ^ ^ | = |n^ 2 ^| = 1. Thus Eqs . 

(4.8.17) and (4.8.18) yield five independent scalar equations as 
follows 


f^(<J>p,0p,<t> 3 ,<J>G'Xf<)>2) — 0 (i — 1,2,. ..,5) (4.8.19) 

Equation system (4.8.19) is the expression in implicit form of 
functions of one variable, for instance of The Theorem of Implicit 

Function System Existence states: 

(i) Consider that a set of parameters 

P = ( <|>p t 9p , ' "X / 4*2 ) (4.8.20) 


satisfies Eq. (4.8.19) and that the Jacobian 

P(f l,f 2 f f 3 f f 4 ,f 5 ) 

D ( <t> p t 9p , <J>g , X , 0 2 ) 


(4.8.21) 


(the partial derivatives are taken at point P°). 

(ii) Then, equation system (4.8.19) can be solved in the 
neighborhood of point P° by functions 


<}> 2 ( $ 1 ) ' <(>p ( <f»i ) / 8 p ( $ ]_ ) , <|>g ( 4 1 1 ) / X ( <|> ^ ) eC (4.8.22) 
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Tne designation means that the functions have continuous 
derivatives at least of the first order. Knowing function <j> 2 (<j>i)r we 
can determine the function of kinematical errors 


a^ 2 ( < j»i ) = ♦ 2 U 1 ) 



(4.3.23) 


4.9 Modification of Generating Surface E p 
Basic Relations 

The proposed method of crowning provides conjugated gear tooth 
surfaces and a localized bearing contact. Using the TCA program it 
can be proven that the kinematical errors caused by the gear 
misalignment are on a very low level. However, the shape of the 
function of kinematic errors is unfavorable, similar to the shape that 
nas been represented in Fig. 1.2a. This disadvantage can be avoided 
by the modification of generating surface £ p . The authors proposed to 
use a surface of revolution instead of a cone surface. Figure 4.9.1 
shows the axial section of generating surface l p that is an arc of the 
circle of radius p. Controlling the magnitude of p, we can provide the 
favorable shape of the function of kinematical errors (Fig. 1.2b) and 
also control tne level of kinematical errors. The solution to this 
problem is based on the relation between the principal curvatures and 
directions for the contacting surfaces and the derivative 


m; 


d 


21 = d^T" ^ m 2 i < f’i^* This relation (Litvin, 1937) is represented as 


follows 
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11 

b 12 

m 

r— 1 


(4.9 

.1) 

21 

b 22 

b 23 

= 0 

31 

b 32 

b 33 




that the 

tooth 

surfaces of 

a crowned pinion and a gear are 

in 

at 

point 

M that coincides 

with the instantaneous center 

of 


rotation (Fig. 4.8.4). It is assumed that the gears are not 
misaligned. In this case the coefficients of determinant (4.9.1) are 
expressed as follows 


b ll b 12 
b 12 b 22 


4 1} + <1 2) 


k (1) + < (2) 

K I K I 


-( 1 ) ( 12 ) 
Sn • a 

0 ( 1 ) ( 12 ) 
e T .a) 


(1)„„(2) 


. ( 2 ) „ w (1) 


( 1 ) n 2„ 


l 21 v ~ A ~2' * v ~f ~f 
(4.9.2) 


_ d * 2 ^ dm 21 _ d d *2 

rn 21 " ' m 21 " d^Y " dij ( d7Y } 

Equations (4.3.11), (4.6.8), (4.6.9) with 0 p = 0 yield (see Fig. 8.4) 


cos t|l 
n = s i n tp 
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( 12 ) _ ( 1 ) ( 2 ) 
0) - 0) “CO 


( 1 ) N 1 
- “ )(1+ NI> k ~f 


wnere N-|_ and N 2 are the number of teeth of the pinion and gear 
R f = ( r i +r 2 > i f 

where r^ and r 2 are the radii of the pitch circles 


( 1) 
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= r lif 
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_ (1) 

•tr 

“ 03 

(2) 

_ (2) 

tr 

— 0) 


l~f 


V.'*:' = U '‘'x(r^'-R f ) = 


bubstitutiny Eqs . (4.9.3) into Eq. (4.9.2), we obtain 
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- (1) . (2) 
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D 13 = 0 


(4.9.3) 
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(4.9.4) 


b 23 ' -“ <1,(1+ ^» 


( 2 ) ^1 ^2 
b 33 = «' 'V (1+ Tq> sin *c + m 21 008 *P 


Equations (4.9.1) and (4.9.4) yield the following relation between the 
aerivative and the principal curvatures 
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(l +!J i ) 2 - (l+gi) sin^r^-qp + qp] 
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( 1 ) ( 2 ) 


21 


N. 


ri c ° s *oTTf I - |C ix > + 'll ’ 1 


(4.9.5) 


Here (see Eqs. (4.6.31) and (4.6.32)) 


• {1) = K (P) 

II II 


.(2) _ 


N l /N 2 


r 1 sin + c 11 r 2 sin,(j c r i sin ^ c 


(4.9.6) 
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(4.9.7) 


Here (Fig. 4.9.1) 


.(P) = _ 1 


II 


where p is the radius of the arc circle that generates the surface of 


revolution 



There are two important particular cases for the application of 

Eq. (4.S.7) (i) The generating surface is a cone. Considering 

that p “ and the generating surface becomes a cone, we obtain that 

' Ifl 2l = ^ the years are not misaligned. In this case, when 

misalignment exists, kinematic errors will occur as shown in Fig. 1.2a 

which are very undersirable . (ii) The generating surface E p is a 

( P ) i ( P ) i 1 

revolute surface, p is positive, <' < 0 ( he' = — ) and 

! i i ■ 1 1 1 p 

m ' 0 (m^ is taken at point M). Then the function of kinematic 
errors will be similar to that as shown in Fig. 1.2b even when the 
years are not aligned. It will be shown below that by choosing an 
appropriate value for p it is possible to keep the kinematical errors 
of misaligned gears on a very low level. 

I 

Surface of Revolution Equations 

I 

The surface of revolution is generated by an arc of circle of 
radius p. The arc parameters provide a common normal for the cone 
surface and the surface of revolution at their point of tangency M 
i (Fig. 4.9.2a). The circular arc is expressed in the auxiliary 
[ coordinate system S Q as follows 


x 0 = p [cos ( ct+e)-cosa] + UpSina 
= u^sinct - 2psin^(a+4) 


y = p [sin ( a + 6) - sina] + 


COSi|> 


cosa 


b o o 

— cosa + 2psin^ cos(ahr) 

cosip 2 2 


z e = 0 


(4.9.9) 
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Tne unit normal to the arc at point M is also the normal to the cone 
surface and surface revolution at this point and it can be expressed 
as follows 


[n ] 
e 


cos( a+6) 
sin( a+6) 
0 


(4.9.10) 


The surface of revolution is generated by rotation of the 
circular arc about the y e -axis and may be represented in coordinate 
system S a as follows 


= t L ael^ e ] 


' [n a J = [Lae] [n Q ] 


Here (Fig. 4.9.2b) 


IL 


ae‘ 


cos 9 p 0 
0 1 
-sine p 0 


sine p 

0 

COS 0 p 


Then we obtain 


X 

a 


- * 6 8 

[UpSina - 2psin-|-sin ( a+-j) ] cos e p 

Y 

a 


cos* COSa + 2psm|cos(a+®) 

Z 

a - 


* 8 8 
[UpSina - 2psin-^sin ( a+^) ] sin 9 p 


in 


cos e p cos (a+6) 
sin( a+6) 
sine cos (a+6) 


(4.9.11) 


(4.9.12) 


(4.9.13) 


(4.9.14) 


53 



Tne installment of tne generating surface in coordinate system S c 
is the same as it was described in Section 4.3. The generating 
surface and its unit normal are represented in S c as follows 

[r< P) ] = [M ca ] [r a ] = (4.9.15) 

» o ft ft - 

2psin [sin( a-^ c ) cos ( a+ -^) - cos 0 p cos ( a-^ c ) sin ( a+ -^) ] 

★ 

+ UpSinacos ( a-^ c ) (cos e p -l ) 

6 8 8 
2psin j [cos ( a-^ c )cos ( a+ -|) + cos 9 p sin( a-t|> ) sin ( a+ -j) ] 

* 

- U p sinasin(a-^ c ) (cos9 p -l) 

* 8 8 
[UpSinct - 2psin sin(a+ Tpisingp 


‘"c^l = 1 L ca' ’ 


(4.9.16) 


cos 9 p cos ( a+ 0)cos ( ) + sin ( <x + 0 ) sin ( a- ip c ) 
-cos9pcos(<x + 9)sin(oi-^ c ) + sin ( a + 0 ) cos ( a-ifi ) 


sin9pCOS ( o+0) 


Pinion Surface 

The pinion surface is represented in coordinate system as 


follows 


f ( 6 1 # 0 p f * p ) 


0 


(4.9.17) 


Lr l J = [M lf J [M fc ] [r c ] 


Here f(8»0p»$p) = 0 is the equation of meshing (see Section 4.4) 
represented as follows 


* 

f(6/9p<f>p) = r 1 sinecos 0 p - u p sinasin ( a + 3 ) ( cos 9 - 1 ) 

- r l( |) p sin(o + 8 )cos(a-i|) ) - cos 9 p cos ( ct+ 3 ) sin ( a-i(» ) = 0 

(4.9.18) 

TCA Program 

Following the procedure of derivations described in Section 4.8, 
we can now develop the TCA program for the case where the generating 
surface £ p is a surface of revolution. Using this program, we can 
determine the kinematical errors caused by the gear. 

Example 


P n =10 N 1 = 20, N 2 =40 C = 1 . 0 1 ( r x +r 2 ) Ay = 5 ' 

$ =20° a = 80° p = 500 in d = 0.176 in 

The results of computation of kinematical errors are shown in 
Fig. 4.9.3. It was found that the function of kinematical errors has 
the desired shape, and the maximal value is 0.16 arc sec. 
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Conclusion 


The authors have developed: (i) The basic principles for the 
optimal geometry of crowned spur gears; (ii) New methods for 
generation of crowned pinion tooth surface; (iii) Tooth contact 
analysis programs for the simulation of meshing and bearing contact. 

The features of the optimal geometry of the crowned spur gears 
are as follows: (i) the gear is provided with a regular involute 
surface; (ii) the pinion tooth surface is deviated from a regular 
involute surface to provide a localized bearing contact and reduce the 
sensitivity of the gears to their misalignment; (iii) the function of 
transmission errors of the misaligned gears is of a parabolic type 
with the small level of maximal errors, less than 5 arc seconds. 

The new methods for generation of the crowned pinion surface need 
simple tools that have to be provided with a generating plane or a 
surface of revolution that slightly deviates from a cone surface. The 
developed TCA programs provide the information on the bearing contact 
and the transmission errors for the misaligned gears. 
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******************************************************* 
* * 


* * 

* PROGRAM I * 

* KINEMATIC ERROR OF A REGULAR INVOLUTE GEAR MESHING * 

* WITH A PINION CROWNED BASED ON PREDESIGNED * 

* KINEMATIC ERROR FUNCTION * 


* * 

* AUTHORS: FAYDOR LITVIN * 

* JIAO ZHANG * 

* * 


******************************************************* 


C 

C PURPOSE 
C 

C THIS PROGRAM IS USED TO CALCULATE THE KINEMATIC ERROR OF A CROWNED 
C PINION MESHING WITH A REGULAR INVOLUTE GEAR. THE PINION IS 

C CROWNED BASED ON PREDESIGNED KINEMATIC ERROR FUNCTION. 

C 

C THIS PROGRAM IS WRITTEN IN FORTRAN77 . IT CAN BEE COMPILED BY V 
C COMPILER IN IBM MAINFRAME OR FORTRAN COMPILER IN VAX SYSTEM. 

C 

IMPLICIT REAL*8(A-H,0-Z) 

DIMENSION Z (99) , ERR (99) , ERROR (99) 

COMMON /BLOCAl/ X (1 1) , Y (10) , A(10, 10) , Y1 (10) , IPVT (10) , WORK (10) , 

+ EPSI, DELTA, NC,NE,NDIM 

COMMON /BL0CA2/ ALPHA, RKSI , RP , RG, RMPG, RL, SA, CA, DR, CK, SK, HDC , C , 

+ S (3, 3) , SF,CF, SAK, CAK, Al ,XP , YP , ZP ,XG, YG, R , CONST (9) 


DEFINE PARAMETERS USED BY PROGRAMS 


(1) IN ANG LP ARE UNIT NUMBERS ASSIGNED TO THE INPUT AND OUTPUT 
DEVICES 

IN-5 
LP-6 

C (2) NDBUG IS USE TO CONTROLL THE AUXILIARY OUTPUT FOR DEBUGGING 
NDBUG-2 

C (3) NC IS THE UPPER LIMITATION OF REPEATATION FOR SOLVING NONLINEAR 
C EQUTIONS ; 

C EPSI IS THE CLEARANCE OF FUNCTION VALUES WHEN THE FUNCTIONS 

C IS CONSIDERED AS SOLVED (ALL FUNTIONS HAVE FORMS OF F (X) -0) ; 

C DELTA IS THE RELATIVE DIFFERENCE FOR TAKING DERIVATIVES 

C NC, EPSI AND DELTA MAY BE CHANGED WHEN SOLUTIONS ARE DIVERGENT 

C OR LESS ACCURATE 

NC-100 
DELTA- l.D-3 
EPSI-l.D-12 

C (4) OTHER PARAMETERS (DON ' T CHANGE) 

NDIM-10 

NE-5 

DR -DAT AN ( 1 . DO) /45 . DO 
C 

C DEFINE INPUT PARAMTERS OF PROBLEM (USE INCH AS UNIT OF LENGTH) 

C 

C (1) PINION AND GEAR: PN-DIAMETRAL PITCH; Nl-PINION TOOTH NUMBER; 

C RMPG- TOOTH NUMBER RATIO (GEAR TOOTH NO./N1); 

C RKSI-PRESSURE ANGLE; HDC-HEIGHT OF DEDENDUM OF PINION; 

C COE-COEFF. OF CENTRAL DISTANCE (USUALLY COE-1.) 
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C R=RADIUS OF ROTATION AXIS FOR GENERATING PINION SURFACE 

PN-10.D0 
Nl=20 
RMPG-2.D0 
RKSI-20.D0*DR 
HDC-1.D0/PN 
COE-l.OlODO 
R-3.5D2 

(2) PRE-DESIGNED FUNCTION: INPUT CONSTANTS HERE AND INPUT FUNCTION 
IN SUBROUTINE FUNC 

MATRIX CONST (J) IS USED FOR CONSTANTS OF FUNCTION 
CONST (1) -2. D0*3. 92975 168 ID-4 

(3) MISALIGNMENT: NMIS-ID NO. (1-CROSSING AXES, 2-INTERSECTING AXES); 
NG-NO. OF MISALIGNED ANGLES TO BE SIMULATED (FROM 0. TO 
(NG-1)*GAMMAI) ; GAMMAI- INCREMENT OF MISALIGNED ANGLE (MINUTE) ; 

NMIS-2 
NG-2 

GAMMAI-5 . DO 

C (4) OUTPUT: FEE 1= INCREMENT OF ROTATION ANGLE OF PINION (DEGREEE) 
FEEI-1 . 0D0*DR 
C 

C DESCRIPTION OF OUTPUT PARAMETERS 
C 

C FEEl-ROTATION ANGLE OF PINION 

C FEE2-ROTATION ANGLE OF GEAR 

C RP-RADIUS OF PINION CONTACT POINT 

C RG= RADIUS OF GEAR CONTACT POINT 

C 

C FIND AUXILIARY VALUES FOR CALCULATION 
RP= FLOAT (Nl) /2 . /PN 
RG=RP*RMPG 
C= (RP+RG) *COE 
CK-DCOS (RKSI) 

SK-DSIN (RKS I) 

NCOEF-360 . DO*DR/FEEI/FLOAT (Nl ) +0 . 3D0 

N=NCOEF*2+l 

NT-N-NCOEF 

FI 1=360. DO/RMPG/FLOAT(N1) 

C DEFINE MATRIX DESCIBEING MISALIGNMENT AND OUTPUT ASSEMBLING CONDITION 
DO 5 1-1,3 
DO 5 J=1 , 3 
S (I , J)=0.D0 

IF (I.EQ.J) S(I,J)=1.DG 
5 CONTINUE 
DO 15 LL-l.NG 

GAMMA=GAMMAI*FLOAT(LL-l) /60.D0*DR 
CG-DCOS (GAMMA) 

SG-DS IN (GAMMA) 

GAMMA= GAMMA / DR * 6 0 . DO 
IF (NMIS.EQ.l) THEN 
WRITE (LP , 10) COE, GAMMA 

10 FORMAT (1H1,///,1X, ' C“ ' ,F4.2, ' * (RP+RG) ; CROSSING ANGLE- ' , 

+ F5 . 1 , ' (M) ' ) 

S(l, 1)-CG 
S (1 , 3) — SG 
S (3 , 1) =SG 
S (3, 3) =CG 
ELSE 
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WRITE (LP, 20) COE, GAMMA 

20 FORMAT (1H1 , /// , IX, ' C- ' , F4.2, ' * (RP+RG) ; INTERSECTING ANGLE- ' , 
+ F5 . 1 , ’ (M) ' ) 

S (2 , 2) -CG 
S (2,3)— SG 
S (3, 2) -SG 
S(3,3)=CG 
END IF 

C SIMULATING MESHING OF PINION AND GEAR (L=l FOR FINDING INITIAL 
C ROTATION CORRESPONDING TO 0 PINION ROTATION) AND OUTPUT RESULTS 
DO 25 L-1,2 
DO 35 1=1,4 
35 X(I)=0.D0 
DO 45 1=1, N 

X (7) =FEEI*FLOAT (I- (N+l) /2) 

IF (L.EQ.l) X (7) =0. DO 
X (5) =DARSIN (RP*X (7) *SK/R) 

X (2) =X (7) /RMPG 
SF=DSIN (X (7) ) 

CF-DCOS (X (7) ) 

CALL NONLIN 
X (8) =X (2) +X (3) 

IF (L.EQ.l) THEN 
XIN=X (8) 

WRITE (LP , 30) 

30 FORMAT (////8X, 'FEEl(D) ' ,8X, ’FEE2(D) ',8X, 'K-ERROR(S) ' ,5X, 

+ 'RP' , 13X, ' RG ' /) 

GO TO 25 
END IF 

X (8) =X (8) -XIN 
X (7) =X (7) /DR 
X (8) =X (8) /DR 

X (9) = (X (8) -X (7) /RMPG) *3600 . DO 
Z (I) =X (8) 

ERR (i) =X (9) 

WRITE (LP , 40) (X(J) , J-7, 11) 

40 FORMAT (IX, 5F15 . 7 , F15 . 7) 

45 CONTINUE 

WRITE (LP , 50) 

50 FORMAT (//,' FIND THE WORKING RANGE FOR ONE TOOTH:'/) 

DO 55 1=1, NT 

X (7) =FEEI*FLOAT (I- (N+l) /2) /DR/2 . DO 

X (8) =X (7) +FII 

KK-I+NCOEF 

ANGLE=Z(KK)-Z(I) 

ERROR (I)-(ANGLE-FII) *3600. DO 
WRITE (LP , 60) X(7) ,X(8) .ANGLE 

60 FORMAT (IX, ’ (' ,F7.2, ’ ’ , F7 . 2 , 1 ) : ’ , F15 . 7 , F15 . 7) 

55 CONTINUE 
DO 65 1=1, NT 
ATEMP2-ERROR (i) 

IF (I.NE.l) THEN 

IF (ATEMP 1 . GT . 0 . DO . AND . ATEMP2 . LE . 0 . DO) GOTO 75 
END IF 

ATEMP1=ATEMP2 
65 CONTINUE 

WRITE (LP , 70) 

70 FORMAT (//IX, ’MESHING IS DISCONTINOUS 1 ) 
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GOTO 25 

75 IF (DABS (ATEMP1) .LT. DABS (ATEMP2) ) 1=1-1 
EMAX-O.DO 
EMIN-0. DO 
DO 85 J-l , NCOEF 
KS-I+J-1 
ET-ERR(KS) 

IF (ET. LT. EMIN) EMIN-ET 
IF (ET.GT.EMAX) EMAX-ET 
85 CONTINUE 

ET-EMAX-EMIN 

KK-I+NCOEF 

WRITE (LP , 80) Z(I) ,Z(KK) ,ET 

80 FORMAT (//IX, 'WORKING RANGE FOR THE GEAR TOOTH: ',F7.2,' ' , 

+ F7.2/1X, 'THE MAXIMUM KINEMATIC ERROR: ',F15.7,' (S)') 

25 CONTINUE 
15 CONTINUE 

WRITE (LP , 90) 

90 FORMAT (1H1) 

STOP 

END 

C 

SUBROUTINE FUNC 
C 

C THIS SUBROUTINE IS USED TO GENERATE FIVE FUNCTIONS, THAT IS, AT 
C CONTACT POINT, POSITION VECTOR AND NORMAL OF PINION AND GEAR 
C MUST COINCIDE 

IMPLICIT REAL*8 (A-H.O-Z) 

COMMON /BLOCAl/ X ( 1 1) , Y (10) , A ( 10 , 10) , Yl (10) , IPVT (10) , WORK (10) , 

+ EPSI, DELTA, NC,NE,NDIM 

COMMON /BLOCA2/ ALPHA, RKSI ,RP , RG, RMPG.RL, SA,CA, DR ,CK, SK,HDC,C, 

+ S (3 , 3) , SF , CF , SAK, CAK, Al ,XP, YP ,ZP,XG, YG, R, CONST (9) 

INPUT PRE-DESIGN ERROR FUNCTION EF AND ITS DIFFEVERTIVE DEF, 
RECALLING X(5) IS FEEP 
EF--CONST (1) *X (5) *X (5) /2.D0 
DEF=-CONST(l)*X(5) 

CT-DCOS (X (A) ) 

ST-DSIN (X (A) ) 

CFEE-DCOS (X (3) ) 

SFEE-DSIN (X (3) ) 

CFEK-DCOS (X(3) +RKSI) 

SFEK-DSIN (X (3) +RKSI) 

GF=DARCOS ( ( 1 . DO+DEF*RMPG/ ( 1 . DO+RMPG) ) *CK) 

CFGK-DCOS (X (5) +GF-RKSI) 

SFGK-DS IN (X (5) +GF-RKS I) 

CFP-DCOS (X (5) ) 

SFP-DSIN (X (5) ) 

CFPG-DCOS (X (5) +GF) 

SFPG-DSIN (X (5) +GF) 

T3=X (5) /RMPG+EF+RKSI-GF 
XI 1=- (RP+RG) *SFP+RG*SFGK+RG*T3*CK’’CFPG 
Yl 1= (RP+RG) *CFP-RG*CFGK+RG*T3*CK*SFPG 
XP-X1 1*CT+R* (CT-1 . DO) 

YP-Y11 

ZP=X11*ST+R*ST 
RNXP=CFPG*CT 
RNYP-SFPG 


68 



no n noon on 


RNZP-CFPG*ST 

XG-RG*SFEE+RG*X(2) *CK*CFEK 
YG=-RG*CFEE+RG*X (2) *CK*SFEK 
ZG=X (1) 

RNXG-CFEK 

RNYG-SFEK 

RNZG-O.DO 

Y (1) »CF*XP+SF*YP-S (1 , 1) *XG-S (1 , 2) *YG~S (1 , 3) *ZG 
Y (2) =CF*YP-SF*XP-S (2,1) *XG~S (2 , 2) *YG-S (2,3) *ZG-C 
Y (3) =ZP-S (3,1) *XG-S (3 , 2) *YG-S (3 , 3) *ZG 

Y (4) «CF*RNYP-SF*RNXP-S (2,1) *RNXG-S (2 , 2) *RNYG-S (2 , 3) *RNZG 
Y (5) =RNZP-S (3,1) *RNXG-S (3 , 2) *RNYG-S (3 , 3) *RNZG 
X (1 0) -DSQRT (XP*XP+ YP*YP) 

X (1 1) =DSQRT (XG*XG+YG*YG) 

WRITE (6,20) XP,YP,ZP,XG,YG,X(6) ,Al,A2,A3 
20 FORMAT (IX, '$$$$' .6F15.7) 

RETURN 

END 

SUBROUTINE NONLIN 

THIS SUBROUTINE IS USED TO SOLVE NONLINEAR EQUATIONS BY NEWTON- 
RAPHSON METHOD 
IMPLICIT REAL*8(A-H,0-Z) 

COMMON /BLOCAl/ X (1 1) , Y (10) , A (10 , 10) , Y1 (10) , IPVT (10) , WORK (10) , 
+ EP S I , DELTA , NC , NE , NDIM 


DO 5 I-l.NC 
CALL FUNC 

WRITE (6,10) I, (X(J) ,Y(J) , J-1,5) 

10 FORMAT (IX, ' *** ' ,15/5 (IX, 2D15 . 7/) ) 

DO 15 J-l.NE 

IF (DABS (Y ( J) ) . GT. EPSI) GO TO 25 
15 CONTINUE 
GO TO 105 
25 DO 35 J=1,NE 
35 Y1(J)-Y(J) 

DO 45 J-l.NE 
DIFF=DABS (X ( J) ) *DELTA 
IF (X(J).EQ.O.DO) DIFF=DELTA 
XMAM-X(J) 

X(J)-X(J)-DIFF 
CALL FUNC 
X (J) =XMAM 
DO 55 K-l.NE 

55 A(K,J)=(Y1(K)-Y(K))/DIFF 
45 CONTINUE 

DO 65 J=1,NE 
65 Y(J)— Y1(J) 

CALL DECOMP (NDIM, NE, A, COND, IPVT, WORK) 
CALL SOLVE (NDIM. NE, A, Y , IPVT) 

DO 75 J=1 ,NE 
75 X(J)=X(J)+Y(J) 

5 CONTINUE 
105 RETURN 
END 
C 

c 
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c 

SUBROUTINE DECOMP (NDIM, N, A, COND, IPVT, WORK) 

C 

IMPLICIT REAL*8(A-H,0-Z) 

DIMENS ION A (NDIM , N) , WORK (N) , IPVT (N) 

C 

C DECOMPOSES AREAL MATRIX BY GAUSSIAN ELIMINATION, 

C AND ESTIMATES THE CONDITION OF THE MATRIX. 

C 

C -COMPUTER METHODS FOR MATHEMATICAL COMPUTATIONS-, BY G. E. FORSYTHE, 
C M. A. MALCOLM, AND C. B. MOLER (PRENTICE-HALL, 1977) 

C 

C USE SUBROUTINE SOLVE TO COMPUTE SOLUTIONS TO LINEAR SYSTEM. 

C 

C IN1UT. . 

C 

C NDIM - DECLARED ROW DIMENSION OF THE ARRAY CONTAINING A 

C N ORDER OF THE MATRIX 

C A MATRIX TO BE TRIANGULAR IZED 

C 

C OUTPUT.. 

C 

C A CONTAINS AN UPPER TRIANGULAR MATRIX U AND A PREMUTED 
C VERSION OF A LOWER TRIANGULAR MATRIX I~L SO THAT 

C (PERMUTATION MATRIX) *A-L*U 

C 

C COND * AN ESTIMATE OF THE CONDITION OF A. 

C FOR THE LINEAR SYSTEM A*X - B , CHANGES IN A AND B 

C MAY CAUSE CHANGES COND TIMES AS LARGE IN X. 

C IF COND+l.O .EQ. COND , A IS SINGULAR TO WORKING 

C PRECISION. COND IS SET TO l.OD+32 IF EXACT 

C SINGULARITY IS DETECTED. 

C 

C IPVT = THE PIVOT VECTOR 

C IPVT (K) = THE INDEX OF THE K-TH PIVOT ROW 

C IPVT (N) - (-1)** (NUMBER OF INTERCHANGES) 

C 

C WORK SPACE.. THE VECTOR WORK MUST BE DECLARED AND INCLUDED 
C IN THE CALL. ITS INlUT CONTENTS ARE IGNORED. 

C ITS OUTPUT CONTENTS ARE USUALLY UNIMPORTANT. 

C 

C THE DETERMINANT OF A CAN BE OBTAINED ON OUTPUT BY 
C DET(A) - IPVT (N) * A(l,l) * A(2,2) * ... * A(N,N) . 

C 

IPVT (N) “1 

IF (N.EQ.l) GO TO 150 
NM1=N-1 

C COMPUTE THE 1-NORM OF A . 

ANORM-O.DO 
DO 20 J-1,N 
T*=0. DO 
DO 10 I-1,N 
10 T=T+DABS (a(I , J) ) 

IF (T.GT.ANORM) ANORM'T 
20 CONTINUE 

C DO GAUSSIAN ELIMINATION WITH PARTIAL 

C PIVOTING. 

DO 70 K-l.NMl 
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KPl-K+l 

C FIND THE PIVOT. 

M-K 

DO 30 I-KPl.N 

IF (DABS (A(I ,K) ) . GT.DABS (A(M,K) ) ) M-I 
30 CONTINUE 
IPVT (K) -M 

IF (M.NE.K) IPVT (N)— IPVT (N) 

T=A (M,K) 

A(M,K)-A(K,K) 

A(K,K)«T 

C SAIP THE ELIMINATION STEP IF PIVOT IS ZERO. 

IF (T.EQ.O.DO) GO TO 70 

COMPUTE THE MULTIPLIERS. 

DO 40 I-KPl.N 
40 A(I,K) — A(I,K)/T 

C INTERCHANGE AND ELIMINATE BY COLUMNS. 

DO 60 J-KPl.N 
T-A(M.J) 

A(M, J)-A(K, J) 

A(K, J)-T 

IF (T.EQ.O.DO) GO TO 60 
DO 50 I-KP1 , N 

50 A ( I , J) “A ( I , J) +A ( I , K) *T 

60 CONTINUE 
70 CONTINUE 
C 

C COND = (1-NORM OF A)*(AN ESTIMATE OF THE 1-NORM OF A-INVERSE) 

C THE ESTIMATE IS OBTAINED BY ONE STEP OF INVERSE ITERATION FOR THE 
C SMALL SINGULAR VECTOR. THIS INVOLVES SOLVING TWO SYSTEMS 
C OF EQUATIONS, (A- TRANSPOSE) *Y - E AND A*Z = Y WHERE E 
C IS A VECTOR OF +1 OR -1 COMPONENTS CHOSEN TO CAUSS GROWTH IN Y. 

C ESTIMATE - (1-NORM OF Z) / (1-NORM OF Y) 

C 

C SOLVE (A-TRANSPOSE) *Y - E . 

DO 100 K«1,N 
T-O.DO 

IF (K.EQ.l) GO TO 90 
KM1=K-1 
DO 80 I-l.KMl 
80 T-T+A(I,K)*WORK(l) 

90 EK” 1 . DO 

IF (T.LT.O.DO) EK— 1.D0 
IF (A(K,K) .EQ.O.DO) GO TO 160 
100 WORK(K) — (EK+T)/A(K,K) 

DO 120 KB“1 ,NMl 
K-N-KB 
T«0.D0 
KPl-K+l 

DO 110 I«KP1,N 
110 T-T+A(I,K)*WORK(K) 

WORK (K) “T 
M-IPVT(K) 

IF (M.EQ.K) GO TO 120 
T-WORK (M) 

WORK (M) -WORK (K) 

WORK (K) -T 
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120 CONTINUE 


C 

YNORM-O . DO 
DO 130 1-1, N 

130 YNORM=YNORM+DABS(WORK(l)) 

C 

C SOLVE A*Z - Y 

CALL SOLVE (NDIM, N, A, WORK, IPVT) 

C 

ZNORM-O.DO 
DO 140 1=1, N 

140 ZNORM=ZNORM+DABS(WORK(l)) 

C 

C ESTIMATE THE CONDITION. 

COND=ANORM*ZNORM/YNORM 
IF (COND.LT. 1 .DO) COND=1.DO 
RETURN 

C 1-BY-l CASE.. 

150 COND= 1 . DO 

IF (A(l, 1) .NE.O.DO) RETURN 
C 

C EXACT SINGULARITY 

160 COND-1.0D32 
RETURN 
END 

SUBROUTINE SOLVE (NDIM, N , A, B , IPVT) 

C 

IMPLICIT REAL*8(A-H,0-Z) 

DIMENSION A (NDIM , N) , B (N) , IPVT (N) 

C 

C SOLVES A LINEAR SYSTEM, A*X = B 

C DO NOT SOLVE THE SYSTEM IF DECOMP HAS DETECTED SINGULARITY. 

C 

C -COMPUTER METHODS FOR MATHEMATICAL COMPUTATIONS-, BY G. E. FORSYTHE, 
C M. A. MALCOLM, AND C. B. MOLER (PRENTICE-HALL, 1977) 

C 

C INlUT. . 

C 

C NDIM = DECLARED ROW DIMENSION OF ARRAY CONTAINING A 
C N ORDER OF MATRIX 

C A TR I ANGULAR IZED MATRIX OBTAINED FROM SUBROUTINE DECOMP 

C B = RIGHT HAND SIDE VECTOR 

C IPVT = PIVOT VECTOR OBTAINED FROM DECOMP 

C 

C OUTPUT.. 

C 

C B = SOLUTION VECTOR, X 
C 

C DO THE FORWARD ELIMINATION. 

IF (N.EQ.l) GO TO 50 
NM1-N-1 
DO 20 K=1 ,NM1 
KP 1 =K+ 1 
M-IPVT(K) 

T=B (M) 

B(M)=B(K) 

B(K)=T 

DO 10 I=KP1,N 
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10 B(I)«B(I)+A(I,K)*T 

20 CONTINUE 

C NOW DO THE BACA SUBSTITUTION. 

DO 40 KB-l.NMl 
KM1-N-KB 
K-KM1+1 

B(K)»B (K) /A(K t K) 

T— B (K) 

DO 30 I-l.KMl 
30 B(I)“B(I)+A(I,K)*T 

40 CONTINUE 
50 B(1)**B(1)/a(1,1) 

RETURN 

END 
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FLOWCHART FOR PROGRAM II 



I 
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ALP-80. D0*DR 
RC-1.D0 

(3) MISALIGNMENT: NMIS-ID NO. (1-CROSSING AXES, 2- INTERSECTING AXES); 
NG-NO. OF MISALIGNED ANGLES TO BE SIMULATED (FROM 0. TO 
(NG-1)*GAMMAI) ; GAMMAI- INCREMENT OF MISALIGNED ANGLE (MINUTE) ; 
N-500 
NG-2 

GAMMAI-5 .DO 

(A) OUTPUT: Z I- INCREMENT OF TOOTH LENTGH OF CROSS SECTION; 

NL-NO. OF CROSS SECTIONS WHERE PINION PROFILE IS SIMULATED; 
N-NO. OF POINTS USED TO DEVIDE TOOL PROFILE IN ANY CROSS 
SECTION (LARGE N, MORE POINTS IS GOT FOR PINION TOOTH PROFILE) 
ZI-0. 10D0 
NL-3 
N-500 
C 

C DESCRIPTION OF OUTPUT PARAMETERS 
C 

C Z 1-DISTANCE BETWEEN CROSS SECTION CONSIDERED AND MIDDLE CROSS 
C SECTION 

C NO-OUTPUT NO. 

C Yl-TOOL SURFACE AUXILIARY VARIABLE 

C XP-X COORDINATE OF PINION PROFILE 

C YP-Y COORDINATE OF PINION PROFILE 

C Rl-RADIUS OF PINION PROFILE 

C FEE-CORESPONDING PINION SURFACE PARAMETER 

C XSH-AVERAGE DEVIATION SHIFT OF CROSS SECTION PROFILE FROM PROFILE 
C OF MIDDLE SECTION (INVOLUTE CURVE) 

C XPE-MAXIMUM DEVIATION OF CROSS SECTION PROFILE FROM INVOLUTE CURVE 
C SDX-STANDARD DEVIATION OF CROSS SECTION PROFILE FROM INVOLUTE 

C CURVE 

C 

C FIND AUXILIARY VALUES FOR CALCULATION 
RP-FLOAT (Nl) /2 . DO/PN 
CL-RC/DSIN (ALP) 

D=CL*DCOS (ALP) 

Al-CL-HDC/DCOS (RKS) 

RPU-RP+HAC 
DO 5 1*1 , NL 
Z- FLOAT (1-1) *ZI 
YY-D-Z/DTAN (ALP) 

WRITE (LP.10) Z 

10 FORMAT (1H1/1X, 'Zl-' .F15.7/1X, ' NO 1 , 10X, ' Y1 ' , 13X, ' XP ' , 13X, ' YP 1 , 

+ 13X, ’R1 ’ ,13X, 'FEE') 

C CALCULATE THE PROPILE OF THE SURFACE CUT BY PLANE Z-CONST 
KKK-0 

DO 15 J-l.N 

Yl=YY*FLOAT(J-l) /FLOAT (N-l) 

F-DSQRT ( 1 . - (Z/ (D- Y 1) /DTAN (ALP) ) **2) 

FEE- ( (D-Yl) *F/DCOS (ALP) -Al* (F*DCOS (ALP) *DCOS (ALP) +DSIN (ALP) *DSIN ( 

+ ALP) ) ) /RP/ (DSIN (ALP) *DCOS (ALP-RKS) -F*DC0S (ALP) "DSIN (ALP-RKS) ) 

XP (I , J) = (D-Yl) * (DTAN (ALP) *F*DCOS (ALP-RKS+FEE) -DSIN (ALP-RKS+FEE) ) 

+ +Al*DSIN (FEE-RKS) -RP*FEE*DCOS (FEE) +RP*DSIN (FEE) 

YP (I , J) — (D-Yl) * (DTAN (ALP) *F*DSIN (ALP-RKS+FEE) +DC0S (ALP-RKS+FEE) ) 

+ +Al*DCOS (FEE-RKS) +RP*FEE*DSIN (FEE) +RP*DC0S (FEE) 

Rl-DSQRT (XP (I , J) **2+YP (I , J) **2) 

IF (Rl.GT.RPU) THEN 
JJ-J-1 
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C. . 
C. . 
c. . 
C.. 
C. . 
C. . 
C. . 
C. . 

c.. 
c.. 
c. . 


******************************************************* 


* * 

* * 

* PROGRAM II * 

* SURFACE OF PINION GENERRATED BY CONE CUTTER * 

* * 

* AUTHORS: FAYDOR LITVIN * 

* JIAO ZHANG * 

* * 

* * 


******************************************************* 


c 

C PURPOSE 
C 

C THIS PROGRAM IS USED TO CALCULATE THE SURFACE OF A PINION WHICH IS 
C GENERATED BY CONE CUTTER 
C 

C NOTE 
C 

C THIS PROGRAM IS WRITTEN IN FORTRAN 77. IT CAN BE COMPILED BY V 
C COMPILER IN IBM MAINFRAME OR FORTRAN COMPILER IN VAX SYSTEM. 

C 

IMPLICIT REAL*8(A-H,0-Z) 

DIMENSION NS (15) ,XP (15 , 999) , YP (15 , 999) ,XERR0R(999) 


DEFINE PARAMETERS USED BY PROGRAMS 


(1) IN ANG LP ARE UNIT NUMBERS ASSIGNED TO THE INPUT AND OUTPUT 

DEVICES 

IN-5 

LP-6 

(2) NDBUG IS USE TO CONTROLL THE AUXILIARY OUTPUT FOR DEBUGGING 
NDBUG-2 

(3) NC IS THE UPPER LIMITATION OF REPEATATION FOR SOLVING NONLINEAR 

EQUTIONS; 

EPSI IS THE CLEARANCE OF FUNCTION VALUES WHEN THE FUNCTIONS 
IS CONSIDERED AS SOLVED (ALL FUNTIONS HAVE FORMS OF F(X)-O); 
RKSTA IS THE RELATIVE DIFFERENCE FOR TAKING DERIVATIVES 
NC, EPSI AND RKSTA MAY BE CHANGED WHEN SOLUTIONS ARE DIVERGENT 
OR LESS ACCURATE 
NC-100 
RKSTA- l.D-3 
EPSI-l.D-12 

(A) OTHER PARAMETERS (DON ' T CHANGE) 

DR-DATAN ( 1 . DO) / A5 . DO 

DEFINE INPUT PARAMTERS OF PROBLEM (USE INCH AS UNIT OF LENGTH) 


(1) PINION AND GEAR: PN-DIAMETRAL PITCH; NI-PINION TOOTH NUMBER; 

RKS-PRESSURE ANGLE (DEGREE) ; HDC-HEIGHT OF DEDENDUM OF PINION; 
HAC-HEIRHT OF ADDENDUM OF PINION 
PN-10.D0 
Nl-20 

RKS=20.D0*DR 
HDC-1.D0/PN 
HAC-1.D0/PN 

C (2) TOOL: ALPHA-HALF OF CONE VERTEX ANGLE (DEGREE) ; 

C RC-RADIUS OF BOTTOM CIRCLE OF CONE 
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XP (I , J) -XP (I , JJ) + (XP (I , J) -XP (I , JJ) ) / (R1-R1TEMP) * (RPU-R1TEMP) 

YP (I , J) -YP (I , JJ) + (YP (I , J) -YP (I , J J) ) / (R1-R1TEMP) * (RPU-R1TEMP) 
FEE-FEETEM+ (FEE-FEETEM) / (R1-R1TEMP) * (RPU-R1TEMP) 

Y1-Y1TEM+ (Y1-Y1TEM) / (Rl-RlTEMP) * (RPU-R1TEMP) 

Rl-DSQRT (XP (I , J) **2+YP (I , J) **2) 

KKK-1 
END IF 

WRITE (LP.20) J,Yl,XP(I,J) , YP (I , J) ,R1, FEE 
20 FORMAT (IX, I4,4F15.7,F15.7) 

IF (KKK.GT.O) GO TO 30 
FEETEM-FEE 
Y1TEM-Y1 
R1TEMP-R1 
15 CONTINUE 
30 NS(I)=J-1 

IF (I.NE.l) GO TO 55 
NS1=NS (1) 

GO TO 5 

C PREPARATION OF INTERPLORATION 
55 NS2=NS(l) 

DO 105 L=1 ,NS2 
J=NS2+1-L 

IF (YP(l.NSl) .GT.YP(I, J)) GO TO 110 
105 CONTINUE 
110 NS2= J 
NREC=2 
XERS=0.D0 
DO 115 L=1 , NS2 
DO 125 J'NREC.NSl 
IF (YP (1 , J) . GT . YP (I , L) ) GO TO 120 
125 CONTINUE 
120 NREC=J 
Jl=J-l 

XERROR (L) =XP (I , L) -XP (1 , Jl) - (YP (I , L) -YP (1 , Jl) ) * (XP (1 , J) -XP (1 , Jl) ) 
+ / (YP (1 , J) -YP (1 , Jl) ) 

XERS=XERS+XERROR (L) 

115 CONTINUE 

XERS =- XERS/ FLOAT (NS2) 

XPE=0.D0 

SDX=0.D0 

IF (NDBUG.GT.2) WRITE (LP,80) 

80 FORMAT (1H1 , 13X, ' NO . ' , 4X, ' DI VI AT ION VALUE') 

DO 45 L-1.NS2 

XERROR (L) “XERROR (L) +XERS 

IF (NDBUG.GT.2) WRITE (LP.50) L.XERROR(L) 

50 FORMAT (13X, 13 , 2F15 . 7) 

IF (DABS (XERROR (L) ) .GT.XPE) XPE=DABS (XERROR (L) ) 

SDX= SDX+XERROR (L ) **2 
45 CONTINUE 

SDX=DSQRT (SDX/FLOAT (NS2) ) 

WRITE (LP.40) XERS ,XPE, SDX 

40 FORMAT (///IX, 'XSH= ' ,E15. 7,5X, ' XPE“ 1 , E15 . 7 , 5X, ' SDX* ' , E15 . 7) 

IF (NDBUG.GT.2) WRITE (LP,60) NS1.NS2 
60 FORMAT (//5X, ' NS1= ' , 16 , 6X, ' NS2= ' ,16) 

5 CONTINUE 
STOP 
END 
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FLOWCHART FOR PROGRAM III 
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C. 

C. 

C. 

C. 

C. 

c. 

c. 

c. 

c. 


c. . . 


******************************************************* 


* * 

* * 

* PROGRAM III * 

* UDDERCUTTING CONDITION FOR PINION GENERRATED * 

* BY CONE CUTTER * 

* * 

* AUTHOR: FAYDOR LITVIN * 

* JIAO ZHANG * 

* * 


c . . * * 

Q ********************* * * * * * * * ****** * * * * ********** * * * * * * * 

C 

C PURPOSE 
C 

C THIS PROGRAM IS USED TO FIND THE UDDERCUTTING CONDITIONS FOR A 
C PINION GENERATED BY CONE CUTTER 
C 

C NOTE 
C 

C THIS PROGRAM IS WRITTEN IN FORTRAN 77. IT CAN BE COMPILED BY V 
C COMPILER IN IBM MAINFRAME OR FORTRAN COMPILER IN VAX SYSTEM. 

C 


IMPLICIT REAL*8(A-H,0-Z) 


DEFINE PARAMETERS USED BY PROGRAMS 


(1) IN ANG LP ARE UNIT NUMBERS ASSIGNED TO THE INPUT AND OUTPUT 
DEVICES 

IN-5 
LP-6 

(2) NDBUG IS USE TO CONTROLL THE AUXILIARY OUTPUT FOR DEBUGGING 
NDBUG-2 

(3) OTHER PARAMETERS (DON ' T CHANGE) 

DR-DATAN ( 1 . DO) / 45 . DO 

DEFINE INPUT PARAMTERS OF PROBLEM (USE INCH AS UNIT OF LENGTH) 

(1) PINION AND GEAR: PN-DIAMETRAL PITCH; Nl-PINION TOOTH NUMBER; 

RKS I “PRESSURE ANGLE; HDC-HEIGHT OF DEDENDUM OF PINION; 

TL-TOOTH LENGTH 

PN-10.D0 
Nl-20 

RKS 1=20. D0*DR 
HDC-1 . DO/PN 
TL-6.D0/PN 

(2) TOOL: ALPHA-HALF OF CONE VERTEX ANGLE (DEGREE) ; 

RC-RADIUS OF BOTTOM CIRCLE OF CONE 

ALPHA-80. DO *DR 
RC-1.D0 

(3) PROBLEM: NPROB-ID NO. OF PROBLEM (-1-GIVEN N1 AND HDC, FIND IF 
UDDERCUTTING OCCUR; 0-GIVEN Nl, FIND MAXIMUM HDC WITHOUT UDDER- 
CUTTING; 1-GIVEN HDC, FIND MINIMUM Nl WITHOUT UDDERCUTTING); 
N-NO. OF THETA VALUES USED CALCULATION (BETWEEN THETAS CORRES- 
PONDING TO MIDDLE SECTION AND EDGE SECTION OF PINION) 

NPROB- 1 
N-ll 
C 


79 


n n n n o o 


DESCRIPTION OF OUTPUT 

OUTPUT IS A STATEMENT BASED ON THE PROBLEM WITHOUT ANY LITERAL 
PARAMETER 

FIND AUXILIARY VALUES FOR CALCULATION 
CL-RC/DS IN (ALPHA) 

SA-DS IN (ALPHA) 

CA-DCOS (ALPHA) 

SK-DSIN(RKSI) 

CK=DCOS (RKSI) 

S AK-DS IN (ALPHA-RKS I ) 

CAK'DCOS (ALPHA-RKS I) 

THEMAX=DARSIN (TL/2 . 0/SA) 

IF (NPROB.NE.O) THEN 
Al= (CL-HDC/CK) *CK 
ACL-Al/CL 
END IF 

IF (NPROB.LT.l) THEN 
RP=FLOAT(Nl)/2./PN 
RCL-RP/CL 
END IF 

WRITE (LP , 90) 

90 FORMAT (1H1) 

IF (NPROB) 5,15,25 
C CHECK IF UNDERCUTTING OCCURS 

5 IF (NDBUG.LT.l) WRITE (LP,10) 

10 FORMAT (1H1/3X, 'NO' ,7X, 'THETA' , 12X, 'A' ,13X, ’ B ' , 14X, ' C ' ,10X, 

+ ' B**2-4.*A*C' ,3X, 'U1/(RC/S IN (ALPHA)) ') 

UMIN-5.D0 
DO 45 I-1,N 

THE-THEMAX*DBLE (FLOAT (I- 1) ) /DBLE (FLOAT (N- 1) ) 

ST-DSIN(THE) 

CT-DCOS (THE) 

C0E1»CA*CAK+CT*SA*SAK 

COE2=SA*CAK-CT*CT*CT*CA*SAK 

COE3=SA*CAK-CT*CA*SAK 

C0E4«ST*ST*SA*CA 

A-COEl*COE2 

B— ACL* (C0E2+C0E1*C0E4) -RCL*COE3**3 

C-ACL*ACL*C0E4 

DD“B*B-4 . *A*C 

IF (DD.GE.O.) THEN 

U1CL” (-B+DSQRT (DD) ) /2 . /A 

IF (UMIN.GT.U1CL) UMIN=U1CL 

U2CL- (-B-DSQRT (DD) ) /2 . /A 

IF (NDBUG.LT.l) WRITE (LP,100) I,THE,A,B,C,DD,U1CL 
100 FORMAT (IX, 14, 8F15 . 7) 

ELSE 

IF (NDBUG.LT.l) WRITE (LP,100) I,THE,A,B,C,DD 
END IF 
45 CONTINUE 

IF (UMIN.LT. l.DO) WRITE (LP.110) 

110 FORMAT (///IX, 'UDDERCUTTING WILL OCCUR FOR YOUR DESIGN') 

IF (UMIN.GE. 1 .DO) WRITE (LP.120) 

120 FORMAT (///IX, 'UDDERCUTTING WILL NOT OCCUR FOR YOUR DESIGN') 
GO TO 35 

C DETERMINE THE MAXIMUM ADDENDUM HEIGHT OF RACA CUTTER 
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15 IF (NDBUG.LT. 1) WRITE (LP,20) 

20 FORMAT (1H1/3X, 'NO' ,7X, 'THETA' , 12X, 'A' , 13X, * B * , 14X, ' C * , 10X, 

+ . b **2-4.*A*C' ,3X, 'ALLOWED RATIO OF HDC/(1/PN)') 

UMIN=5.D0 
DO 55 1=1, N 

THE=THEMAX*DBLE (FLOAT (I- 1) ) /DBLE (FLOAT (N-l) ) 

ST-DSIN(THE) 

CT-DCOS (THE) 

C0E1=CA*CAK+CT*SA*SAK 

COE2»SA*CAK-CT*CT*CT*CA*SAK 

C0E3=SA*CAK-CT*CA*SAK 

C0E4=ST*ST*SA*CA 

C“COEl*COE2-RCL*COE3**3 

B=-(COE2+COEl*COE4) 

A=COE4 

DD=B*B-4.*A*C 
IF (DD.GE.O.) THEN 
IF (A.NE.O.) THEN 
U1CL= (-B-DSQRT (DD) ) /2 . /A 
U2CL= (-B+DSQRT (DD) ) /2 . /A 
ELSE 

U1CL=-C/B 
U2CL=0 . 

END IF 

U1CL= (CK-U1CL) *CL*PN 
U2CL= (CK-U2CL) *CL*PN 

IF (NDBUG.LT. 1) WRITE (LP,100) I , THE, A, B , C , DD, U1CL 
ELSE 

IF (NDBUG.LT. 1) WRITE (LP,100) I , THE , A, B , C , DD 
END IF 

IF (UMIN.GT.U1CL) UMIN=U1CL 
55 CONTINUE 

WRITE (LP.200) UMIN 

200 FORMAT (///IX, 'TO AVOID UDDERCUTTING, IT IS NECESSARY TO KEEP DEDE 
+NDUM OF PINION <= ' , F10. 7, '/PN ' ) 

GO TO 35 

C DETERMINE THE MINIMUM NO. OF TEETH FOR UNUNDERCUTTING 
25 IF (NDBUG.LT. 1) WRITE (LP,30) 

30 FORMAT (1H1/3X, 'NO' ,7X, 'THETA' ,7X, 'NO. OF TEETH') 

RNMAX=0. 

DO 65 1=1, N 

THE=THEMAX*DBLE (FLOAT (I- 1) ) /DBLE (FLOAT (N-l) ) 

ST=DSIN (THE) 

CT=DCOS (THE) 

C0E1=CA*CAK+CT*SA*SAK 

COE2=SA*CAK-CT*CT*CT*CA*SAK 

COE3=SA*CAK-CT*CA*SAK 

C0E4=ST*ST*SA*CA 

RR= (COEl*COE2-ACL 5 ' (COE2+COEl*COE4) +ACL*ACL*C0E4) /COE3**3 
RN=2.*RR*PN*CL 
IF (RNMAX.LT. RN) RNMAX=RN 
IF (NDBUG.LT. 1) WRITE (LP,100) I , THE , RN 
65 CONTINUE 

WRITE (LP.300) RNMAX 

300 FORMAT (/// IX, 'WITHOUT UDDERCUTTING, MINIMUM TOOTH NO. OF PINION 
+IS: ' , Fll . 7) 

35 WRITE (LP.400) 

400 FORMAT (1H1) 
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******************************************************* 


PROGRAM IV 

CONTACT ELLIPSIS FOR PINION GENERATED BY CONE 
CUTTER IN MESHING WITH REGULAR GEAR 

AUTHORS: FAYDOR LITVIN 
JIAO ZHANG 


* * * * * * * * ******************* * * * * * * * * * * * * * * * * * * * * * * * * * * * * 


PURPOSE 

THIS PROGRAM IS USED TO CALCULATE THE CONTACT ELLIPSIS WHEN PINION 
GENERATED BY CONE CUTTER IN MESHING WITH REGULAR INVOLUTE GEAR 

NOTE 

THIS PROGRAM IS WRITTEN IN FORTRAN 77. IT CAN BE COMPILED BY V 
COMPILER IN IBM MAINFRAME OR FORTRAN COMPILER IN VAX SYSTEM. 


IMPLICIT REAL*8(A-H,0-Z) 


DEFINE PARAMETERS USED BY PROGRAMS 


(1) IN ANG LP ARE UNIT NUMBERS ASSIGNED TO THE IN1UT AND OUTPUT 

DEVICES 

IN=5 

LP-6 

(2) NDBUG IS USE TO CONTROLL THE AUXILIARY OUTPUT FOR DEBUGGING 
NDBUG=2 

(3) OTHER PARAMETERS (DON’T CHANGE) 

DR-DATAN ( 1 . DO) / 45 . DO 

DEFINE INPUT PARAMTERS OF PROBLEM (USE INCH AS UNIT OF LENGTH) 

(1) PINION AND GEAR: PN-DIAMETRAL PITCH; Nl=PINI0N TOOTH NUMBER; 

RMPG-TOOTH NUMBER RATIO (GEAR TOOTH NO./N1); 

RKSI-PRESSURE ANGLE 
PN-10.D0 
Nl-20 
RMPG=2 . DO 
RKSI=20.D0*DR 

(2) TOOL: ALPHA=HALF OF CONE VERTEX ANGLE (DEGREE) ; 

RC=RADIUS OF BOTTOM CIRCLE OF CONE; 

ALPHA=89 . 5D0*DR 
RC= 1 . DO 

(3) DEFORMATION: DEL=CONTACT DEFORMATION AT CONTACT POINT 
DEL=4.D-4 

(4) OUTPUT: PD=RECIPROCAL OF INCREMENT OF S VALUE (S IS THE DISTANCE 

CONTACT POINT AND INSTANTANEOUS CENTER) 

PD=100.D0 


DESCRIPTION OF OUTPUT PARAMETER 


S=DISTANCE BETWEEN CONTACT POIN" JJD INSTANTANEOUS CENTER 
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C Rl-PINION RADIUS OF CONTACT POINT 

C FEE-PINION ROTATION ANGLE CORRESPONDING TO CONTACT POINT 
C A-LENGTH OF HALF SHORT AXIS OF CONTACT ELLIPSE 
C B-LENTH OF HALF LONG AXIS OF CONTACT ELLIPSE (ALONG DIRECTION OF 
C GEAR TOOTH LEHGTH) 

C 

C FIND AUXILIARY VALUES FOR CALCULATION 
RP-FLOAT (Nl)/2./PN 
SA-DS IN (ALPHA) 

CA-DCOS (ALPHA) 

SK-DSIN(RKSI) 

CK-DCOS(RKSI) 

RL-RC/SA 

SMAX“ (DSQRT ( ( (1 . DO+2 . DO/FLOAT (Nl) ) /CK) **2-1 . DO) -SK/CK) *CK*RP 
NA= IDINT (SHAX*PD+0 . 5D0) 

SMAX«FLOAT(NA)/PD 

SMIN- (DSQRT ( ( (1 . +2 . /FLOAT (Nl) /RMPG) /CK) **2-1 . ) -SK/CK) *CK*RP*RMPG 
NB— IDINT (SMIN*PD+0 . 5D0) 

SMIN“FLOAT(NB)/PD 

N-l+NA-NB 

WRITE (LP , 10) RP,N1,RMPG 

10 FORMAT (1H1.///5X, 'RADIUS OF THE PITCH CIRCLE OF PINION (RP) : ' , 
+F15.7//5X, 'TOOTH NO. OF PINION (Nl) I8//5X, 1 RATIO OF ANGULAR VELO 
+CITY (OMEGA OF PINION/OMEGA OF GEAR) : 1 , F15 . 7///1X, 1 OUTPUT DESIGNAT 
+ ION: S=DISTANCE BETWEEN CONTACT POINT AND INSTANTANEOUS CENTER'/ 
+21X, 1 R1=PINI0N RADIUS OF CONTACT POINT; FEE-PINION ROTATION CORRES 
+PONDING TO CONTACT POINT 1 /2 IX, ' A= LENGTH OF HALF SHORT AXIS; B=LENG 
+TH OF HALF LONG AXIS (ALONG DIRECTION OF GEAR TOOTH LEHGTH) ’///// 
+10X, 'S’.IAX, 'Rl’.lOX, 'FEEl(D) ' , 1 IX, 'A' , 14X, ' B ' , 13X, ' B/A '//) 

C FIND POSITION AND DIMENSION OF CONTACT ELLIPSE 
DO 5 L=1 , N 

S-SMIN+ (SMAX-SMIN) /FLOAT (N-l) *FLOAT (L-l) 

A=DSQRT (DEL*RP* (SK+S/RP) * (RMPG*SK~S/RP) *2 . DO/ (1 . DO+RMPG) /SK) 
B-DSQRT (DEL*2. DO*SA/CA* (RL-2 . DO*RP/CK/FLOAT (Nl) -S*SK/CK) ) 

C=B/A 

FEE-S/RP/CK 

R1-RP*DSQRT (1 . D0+ (FEE*CK) **2+2 . DO*FEE*CK*SK) 

FEE=FEE/DR 

WRITE (LP , 20) S , R1 , FEE , A, B , C 
20 FORMAT (1X.6F15.7) 

5 CONTINUE 
WRITE (LP.30) 

30 FORMAT (1H1) 

STOP 

END 
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c. . 
c. . 
c. . 
c.. 
c. . 
c. . 
c.. 
c.. 
c. . 
c. . 
c.. 
c. . 
c 

C PURPOSE 

THIS PROGRAM IS USED TO CALCULATE THE KINEMATIC ERROR OF A PINION 
GENERATEED BY CONE CUTTER MESHING WITH A MISALIGNED REGULAR 
INVOLUTE GEAR 

C NOTE 
C 

C THIS PROGRAM IS WRITTEN IN FORTRAN 77. IT CAN BE COMPILED BY V 
C COMPILER IN IBM MAINFRAME OR FORTRAN COMPILER IN VAX SYSTEM. 

C 

IMPLICIT REAL*8(A-H,0-Z) 

DIMENSION Z (99) , ERR (99) .ERROR (99) 

COMMON /BLOCAl/ X (1 1) , Y (10) ,A(10, 10) , Y1 (10) , IPVT (10) , WORK (10) , 

+ EPSI, DELTA, NC.NE.NDIM 

COMMON /BLOCA2/ ALPHA, RKSI,RP,RG,RMPG,RL,SA,CA, DR, CK.SK.HDC, 

+ C, S (3 , 3) , SF , CF , SAK, CAK, Al ,XP, YP,ZP ,XG, YG,R 

DEFINE PARAMETERS USED BY PROGRAMS 

(1) IN ANG LP ARE UNIT NUMBERS ASSIGNED TO THE INPUT AND OUTPUT 
DEVICES 

IN=5 
LP-6 

(2) NDBUG IS USE TO CONTROLL THE AUXILIARY OUTPUT FOR DEBUGGING 
NDBUG=2 

(3) NC IS THE UPPER LIMITATION OF REPEATATION FOR SOLVING NONLINEAR 
EQUTIONS ; 

EPSI IS THE CLEARANCE OF FUNCTION VALUES WHEN THE FUNCTIONS 
IS CONSIDERED AS SOLVED (ALL FUNTIONS HAVE FORMS OF F(X)=0); 
DELTA IS- THE RELATIVE DIFFERENCE FOR TAKING DERIVATIVES 
NC, EPSI AND DELTA MAY BE CHANGED WHEN SOLUTIONS ARE DIVERGENT 
OR LESS ACCURATE 
NC=100 
DELTA=l.D-3 
EPSI*=1 .D-12 

(4) OTHER PARAMETERS (DON ' T CHANGE) 

NDIM=10 
NE=5 

DR=DATAN ( 1 . DO) /45 . DO 

DEFINE INPUT PARAMTERS OF PROBLEM (USE INCH AS UNIT OF LENGTH) 

(1) PINION AND GEAR: PN=DIAMETRAL PITCH; Nl-PINION TOOTH NUMBER; 
RMPG- TOOTH NUMBER RATIO (GEAR TOOTH NO./N1); 

RKSI=PRESSURE ANGLE; HDC=HEIGHT OF DEDENDUM OF PINION; 
COE=COEFF. OF CENTRAL DISTANCE (USUALLY COE=l.) 



******************************* ******* ***************** 


* * 

* * 

* PROGRAM V * 

* KINEMATIC ERROR OF THE PINION GENERATED BY CONE * 

* CUTTER MESHING WITH A MISALIGNED REGULUAR GEAR * 

* * 

* AUTHORS: FAYDOR LITVIN * 

* JIAO ZHANG * 

* * 

* * 


******************************************************* 
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PN-10.D0 

Nl-20 

RMPG-2.D0 

RKSI=20.D0*DR 

HDC-l.DO/PN 

COE-1.005D0 

C (2) TOOL: ALPHA-HALF OF CONE VERTEX ANGLE (DEGREE) ; 

C RC-RADIUS OF BOTTOM CIRCLE OF CONE 

C R-RADIUS OF ARC (FOR EXACT CONE CHOOSE R >1 .D6*Nl/2/PN) 

ALPHA-80. D0*DR 

RC-l.DO 

R-1.0D3 

C (3) MISALIGNMENT: NMIS-ID NO . (1-CROSSING AXES , 2=INTERSECTING AXES) ; 

C NG=NO. OF MISALIGNED ANGLES TO BE SIMULATED (FROM 0. TO 

C (NG-1)*GAMMAI) ; GAMMAI= INCREMENT OF MISALIGNED ANGLE (MINUTE) ; 

NMIS=1 
NG-2 

GAMMAI-5.D0 

C (4) OUTPUT: FEE I = INCREMENT OF ROTATION ANGLE OF PINION (DEGREEE) 
FEEI=1.0D0*DR 
C 

C DESCRIPTION OF OUTPUT PARAMERTERS 
C 

C FEE1-ROTATION ANGLE OF PINION 

C FEE2=ROTATION ANGLE OF GEAR 

C RP-RADIUS OF PINION CONTACT POINT 

C RG-RADIUS OF GEAR CONTACT POINT 

C 

C FIND AUXILIARY VALUES FOR CALCULATION 
RP-FLOAT (Nl)/2./PN 
RG-RP*RMPG 
C= (RP+RG) *COE 
CA-DCOS (ALPHA) 

SA=DS IN (ALPHA) 

CK-DCOS(RKSI) 

SK=DSIN (RKSl) 

SAK=DS IN (ALPHA-RKS I ) 

CAK=DCOS (ALPHA-RKS I ) 

RL=RC/SA 

A1=RL-HDC/CK 

NCOEF-360 . D0*DR/FEEI/ FLOAT (Nl) +0. 3D0 

N=NCOEF*2+l 

NT=N-NCOEF 

FII=360.D0/RMPG/FLOAT (Nl) 

C DEFINE MATRIX DESCIBEING MISALIGNMENT AND OUTPUT ASSEMBLING CONDITION 
DO 5 1-1,3 
DO 5 J=l,3 
S(I, J)-0.D0 

IF (I.EQ.J) S (I , J)=l .DO 
5 CONTINUE 

DO 15 LL-l.NG 

GAMMA- GAMMA I * F L 0 AT (LL~0) /60 . D0*DR 
CG-DCOS (GAMMA) 

SG-DS IN (GAMMA) 

GAMMA-GAMMA/DR *60 . DO 
IF (NMIS.EQ. 1) THEN 
WRITE (LP , 10) COE, GAMMA 

10 FORMAT (1H1.///.1X, ' C= ' ,F5.3, '* (RP+RG) ; CROSSING ANGLE- ' , 


c-> 
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+ F5 . 1 , ’ (M) ' ) 

S(1,1)-CG 
S (1 , 3) — SG 
S(3,1)-SG 
S(3,3)-CG 
ELSE 

WRITE (LP,20) COE, GAMMA 

20 FORMAT (1H1 , /// , IX, ' C“ ' ,F4.2, ' * (RP+RG) ; INTERSECTING ANGLE- ' , 
+ F5 . 1 , ' (M) ' ) 

S (2, 2) =CG 
S(2,3)“-SG 
S(3,2)=SG 
S(3,3)-CG 
END IF 

C SIMULATING MESHING OF PINION AND GEAR (L“l FOR FINDING INITIAL 
C ROTATION CORRESPONDING TO 0 PINION ROTATION) AND OUTPUT RESULTS 
DO 25 L=1 , 2 
DO 35 1=1,4 
35 X(I)=0.D0 
DO 45 1=1, N 

X (7) =FEEI*FLOAT (I- (N+l) /2) 

IF (L.EQ.l) X(7)=0.D0 
X (5) -DARSIN (RP*X (7) *SK/R) 

X (2) =X (7) /RMPG 
SF=DSIN (X (7) ) 

CF*DCOS (X (7) ) 

CALL NONLIN 
X(8)=X(2)+X(3) 

IF (L.EQ.l) THEN 
XIN=X (8) 

WRITE. (LP,30) 

30 FORMAT (////8X, ' FEE1 (D) 1 , 8X, ' FEE2 (D) ' ,8X, ' K- ERROR (S) ' , 5X, 

+ ’ RP ' ,13X, 'RG'/) 

GO TO 25 
END IF 

X (8) “X (8) -XIN 
X (7) =X (7) /DR 
X(8)=X(8)/DR 

X (9) = (X (8) -X (7) /RMPG) *3600 . DO 
Z (I) =X (8) 

ERR (I) =X (9) 

WRITE (LP.40) (X(J) , J-7, 11) 

40 FORMAT (IX, 5F15 . 7 , F15 . 7) 

45 CONTINUE 

WRITE (LP , 50) 

50 FORMAT (//,' FIND THE WORKING RANGE FOR ONE TOOTH: 7) 

DO 55 1=1, NT 

X(7)=FEEI*FL0AT(l-(N+l)/2)/DR/2.D0 

X(8)=X(7)+FII 

RK=I+NC0EF 

ANGLE=Z(KK)-Z(I) 

ERROR (I) = (ANGLE-FII) *3600 . DO 
WRITE (LP.60) X(7) ,X(8) .ANGLE 

60 FORMAT (IX, ’ (' ,F7.2, ’ ' ,F7.2, ') : ’ , F15 . 7 , F15 .7) 

55 CONTINUE 
DO 65 1=1, NT 
ATEMP2=ERR0R (I) 

IF (I.NE.l) THEN 
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IF (ATEMP1 . GT. O.DO. AND. ATEMP2 . LE. O.DO) GOTO 75 
END IF 

ATEMP1-ATEMP2 
65 CONTINUE 

WRITE (LP , 70) 

70 FORMAT (//IX, ’MESHING IS DISCONTINOUS ' ) 

GOTO 25 

75 IF (DABS (ATEMP 1) .LT. DABS (ATEMP2)) 1=1-1 
EMAX-O.DO 
EMIN=0. DO 
DO 85 J-l.NCOEF 
KS-I+J-1 
ET-ERR (KS) 

IF (ET. LT.EMIN) EMIN*ET 
IF (ET.GT.EMAX) EMAX-ET 
85 CONTINUE 

ET-EMAX-EMIN 

KK=I+NCOEF 

WRITE (LP , 80) Z (I) , Z (KK) ,ET 

80 FORMAT (//IX, 'WORKING RANGE FOR THE GEAR TOOTH: ’,F7.2,’ 

+ F7. 2/ IX, 'THE MAXIMUM KINEMATIC ERROR: ’ , F15 . 7 , ' (S)') 

25 CONTINUE 
15 CONTINUE 

WRITE (LP,90) 

90 FORMAT (1H1) 

STOP 

END 

SUBROUTINE FUNC 

THIS SUBROUTINE IS USED TO GENERATE FIVE FUNCTIONS, THAT IS, AT 
CONTACT POINT, POSITION VECTOR AND NORMAL OF PINION AND GEAR 
MUST COINCIDE 
IMPLICIT REAL*8(A-H,0-Z) 

COMMON /BLOCAl/ X (1 1) , Y (10) , A (10 , 10) , Y1 (10) , IPVT (10) , WORK (10) , 
+ EPS I , DELTA , NC , NE , NDIM 

COMMON /BLOCA2/ ALPHA, RKSI ,RP , RG, RMPG, RL, SA, CA, DR , CK, SK, HDC , 

+ C, S (3 , 3) , SF.CF, SAK, CAK, Al ,XP , YP , ZP ,XG, YG, R 

CT*DCOS (X (4) ) 

ST*DSIN (X (4) ) 

CFEE=DCOS (X (3) ) 

SFEE”DSIN (X (3) ) 

CFEK“DCOS (X (3) +RKSI) 

SFEK“DSIN (X (3) +RKSI) 

CAB=DCOS (ALPHA+X (5) ) 

SAB=DSIN (ALPHA+X (5)) 

CAB2“DCOS (ALPHA+X (5) /2 . DO) 

SAB2=DS IN (ALPHA+X (5) /2 . DO) 

SB”DSIN (X (5) ) 

SB2=DSIN (X (5) /2 . DO) 

A2=Al*SA*(CT-l) 

X (6) * (A2*SAB-R*SB*CT) / (SAB*CAK-CT*CAB*SAK) /RP 
CFP“DC0S (X (6) ) 

SFP=DSIN (X (6) ) 

CAKF-DCOS (ALPHA-RKSI+X (6) ) 

SAKF-DSIN (ALPHA-RKSI+X (6) ) 

A3-2.D0*R*SB2 

T3-CAB2*SAKF-CT*SAB2*CAKF 
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XP-A3* (CAB2*SAKF-CT*SAB2*CAKF) +A2*CAKF-RP*X (6) *CFP+RP*SFP 

YP-A3* (CAB2*CAKF+CT*SAB2*SAKF) -A2*SAKF+RP*X (6) *SFP+RP*CFP 

ZP“ (Al*SA-A3*SAB2) *ST 

RNXP-CT*CAB*CAKF+SAB*SAKF 

RKYP— CT*CAB*SAKF+SAB*CAKF 

RNZP-ST*CAB 

XG“RG*SFEE+RG*X (2) *CK*CFEK 
YG—RG*CFEE+RG*X (2) *CK*SFEK 
ZG“X (1) 

RNXG-CFEK 

RNYG-SFEK 

RNZG-O.DO 

T1-CF*XP+SF*YP 

T2“S (1 , 1) *XG+S (1 , 2) *YG+S (1 , 3) *ZG 
WRITE (6,50) T1,T2,XP,YP,X(6) 

50 FORMAT (IX, ' Tl= ' , F15 . 7 , 5X, ' T2= ' , F15 . 7 , 2F15 . 7) 

Y (1)=CF*XP+SF*YP-S (1 , 1) *XG-S (1 , 2) *YG-S (1 , 3) *ZG 
Y (2) =CF*YP-SF*XP-S (2,1) *XG-S (2 , 2) *YG~S (2 , 3) *ZG-C 
Y (3) =ZP-S (3,1) *XG-S (3 , 2) *YG-S (3 , 3) *ZG 

Y (4) =CF*RNYP-SF*RNXP-S (2 , 1) *RNXG~S (2 , 2) *RNYG~S (2 , 3) *RNZG 
Y (5) =RNZP-S (3,1) *RNXG-S (3 , 2) *RNYG-S (3 , 3) *RNZG 
X (10) -DSQRT (XP*XP+YP*YP) 

X (1 1) “DSQRT (XG*XG+YG*YG) 

WRITE (6,20) XP,YP,ZP,XG,YG,X(6) ,A1,A2,A3 
20 FORMAT (IX, ’$$$$' .6F15.7) 

RETURN 

END 

SUBROUTINE NONLIN 

THIS SUBROUTINE IS USED TO SOLVE NONLINEAR EQUATIONS BY NEWTON- 
RAPHSON METHOD 
IMPLICIT REAL*8(A-H,0-Z) 

COMMON /BLOCAl/ X (11) , Y (10) ,A(10, 10) , Y1 (10) , IPVT(10) , WORK(IO) , 
+ EPSI, DELTA, NC,NE,NDIM 


DO 5 I-1,NC 
CALL FUNC 

WRITE (6,10) I, (X(J) ,Y(J) ,J=1,5) 
10 FORMAT (IX, '***' ,I5/5(1X,2D15.7/)) 
DO 15 J-1,NE 

IF (DABS (Y ( J) ) .GT. EPSI) GO TO 25 
15 CONTINUE 
GO TO 105 
25 DO 35 J-l.NE 
35 Y1(J)-Y(J) 

DO 45 J-l.NE 
DIFF-DABS (X ( J) ) *DELTA 
IF (X(J) .EQ.O.DO) DIFF-DELTA 
XMAM-X(J) 

X (J) -X (J) -DIFF 
CALL FUNC 
X (J) -XMAM 
DO 55 K-l.NE 

55 A(K, J)-(Y1(K)-Y(K))/DIFF 
45 CONTINUE 
DO 65 J-l.NE 
65 Y (J) — Y1 (J) 


91 


non 


CALL DECOMP (NDIM,NE,A,COND, IPVT.WORK) 
CALL SOLVE (NDIM,NE,A, Y, IPVT) 

DO 75 J-l.NE 
75 X(J)-X(J)+Y(J) 

5 CONTINUE 
105 RETURN 
END 


SUBROUTINE DECOMP (NDIM,N,A ( COND, IPVT.WORK) 

C 

IMPLICIT REAL*8(A-H,0-Z) 

DIMENSION A(NDIM,N) , WORK (N) , IPVT (N) 

C 

C DECOMPOSES AREAL MATRIX BY GAUSSIAN ELIMINATION, 

C AND ESTIMATES THE CONDITION OF THE MATRIX. 

C 

C -COMPUTER METHODS FOR MATHEMATICAL COMPUTATIONS-, BY G. E. FORSYTHE, 
C M. A. MALCOLM, AND C. B. MOLER (PRENTICE-HALL, 1977) 

C 

C USE SUBROUTINE SOLVE TO COMPUTE SOLUTIONS TO LINEAR SYSTEM. 

C 

C IN1UT. . 

C 

C NDIM - DECLARED ROW DIMENSION OF THE ARRAY CONTAINING A 

C N ORDER OF THE MATRIX 

C A MATRIX TO BE TRIANGULARIZED 

C 

C OUTPUT.. 

C 

C A CONTAINS AN UPPER TRIANGULAR MATRIX U AND A PREMUTED 
C VERSION OF A LOWER TRIANGULAR MATRIX I-L SO THAT 

C (PERMUTATION MATRIX) *A-L*U 

C 

C COND - AN ESTIMATE OF THE CONDITION OF A. 

C FOR THE LINEAR SYSTEM A*X - B , CHANGES IN A AND B 

C MAY CAUSE CHANGES COND TIMES AS LARGE IN X. 

C IF C0ND+1.0 .EQ. COND , A IS SINGULAR TO WORKING 

C PRECISION. COND IS SET TO 1.0D+32 IF EXACT 

C SINGULARITY IS DETECTED. 

C 

C IPVT - THE PIVOT VECTOR 

C IPVT (K) - THE INDEX OF THE K-TH PIVOT ROW 

C IPVT (N) - (-1)** (NUMBER OF INTERCHANGES) 

C 

C WORK SPACE.. THE VECTOR WORK MUST BE DECLARED AND INCLUDED 
C IN THE CALL. ITS IN1UT CONTENTS ARE IGNORED. 

C ITS OUTPUT CONTENTS ARE USUALLY UNIMPORTANT. 

C 

C THE DETERMINANT OF A CAN BE OBTAINED ON OUTPUT BY 
C DET(A) = IPVT (N) * A (1 , 1) * A(2,2) * ... * A(N,N) . 

C 

IPVT (N) =1 

IF (N.EQ.l) GO TO 150 
NM1-N-1 

C COMPUTE THE 1-NORM OF A . 

ANORM=0 . DO 
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DO 20 J-1,N 
T=0.D0 
DO 10 I-1,N 
10 T«T+DABS(A(I, J)) 

IF (T.GT.ANORM) ANORM-T 
20 CONTINUE 

DO GAUSSIAN ELIMINATION WITH PARTIAL 
PIVOTING. 

DO 70 K-l.NMl 
KP1-K+1 

FIND THE PIVOT. 

M-K 

DO 30 I-KPl.N 

IF (DABS (A(I,K) ) .GT.DABS (A(M,K) )) M-I 
30 CONTINUE 
IPVT(K) -M 

IF (M.NE.K) IPVT(N) — IPVT(N) 

T-A(M,K) 

A(M,K)=A(K,K) 

A(K,K)-T 

SAIP THE ELIMINATION STEP IF PIVOT IS ZERO. 
IF (T.EQ.O.DO) GO TO 70 

COMPUTE THE MULTIPLIERS. 

DO 40 I-KPl.N 
40 A(I,K>— A(I,K)/T 

C INTERCHANGE AND ELIMINATE BY COLUMNS. 

DO 60 J-KPl.N 
T-A(M, J) 

A(M,J)-A(K,J) 

A(K,J)=T 

IF (T.EQ.O.DO) GO TO 60 
DO 50 I-KPl.N 

50 A(I,J)-A(I,J)+A(I,K)*T 

60 CONTINUE 
70 CONTINUE 
C 

C COND - (1-NORM OF A) * (AN ESTIMATE OF THE 1-NORM OF A- INVERSE) 

C THE ESTIMATE IS OBTAINED BY ONE STEP OF INVERSE ITERATION FOR THE 
C SMALL SINGULAR VECTOR. THIS INVOLVES SOLVING TWO SYSTEMS 
C OF EQUATIONS, (A- TRANSPOSE) *Y = E AND A*Z - Y WHERE E 
C IS A VECTOR OF +1 OR -1 COMPONENTS CHOSEN TO CAUSS GROWTH IN Y. 

C ESTIMATE - (1-NORM OF Z)/ (1-NORM OF Y) 

C 

C SOLVE (A- TRANSPOSE) *Y = E . 

DO 100 K=1,N 
T=0.D0 

IF (K.EQ.l) GO TO 90 
KM1-K-1 
DO 80 I-1.KM1 
80 T-T+A(I,K)*WORK(I) 

90 EK-1.D0 

IF (T.LT.O.DO) EK— 1.D0 
IF (A(K,K) .EQ.O.DO) GO TO 160 
100 WORK (K)=-(EK+T) /A(K,K) 

DO 120 KB-l.NMl 
K-N-KB 
T=0.D0 
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KP1-K+1 

DO 110 I-KPl.N 
110 T-T+A(I,K)*WORK(K) 

WORK(K) -T 
M-IPVT(K) 

IF (M.EQ.K) GO TO 120 
T-WORK(M) 

WORK (M) -WORK (K) 

WORK (K) =T 
120 CONTINUE 
C 

Y NORM- 0. DO 
DO 130 1-1 , N 

130 YNORM«YNORM+DABS(WORK(I)) 

c 

C SOLVE A*Z - Y 

CALL SOLVE (NDIM,N,A,WORK, IPVT) 

C 

ZNORM-O. DO 
DO 140 1-1 ,N 

140 ZNORM-ZNORM+DABS (WORK(l) ) 

C 

C ESTIMATE THE CONDITION. 

COND=ANORM*ZNORM/YNORM 
IF (COND.LT. l.DO) COND-l.DO 
RETURN 

C 1-BY-l CASE.. 

150 COND-l.DO 

IF (A ( 1 , 1) .NE.O.DO) RETURN 
C 

C EXACT SINGULARITY 

160 COND-1 . 0D32 
RETURN 
END 

SUBROUTINE SOLVE (NDIM.N, A, B , IPVT) 

C 

IMPLICIT REAL*8(A-H,0-Z) 

DIMENSION A(NDIM,N) ,B (N) , IPVT(N) 

SOLVES A LINEAR SYSTEM, A*X - B 

DO NOT SOLVE THE SYSTEM IF DECOMP HAS DETECTED SINGULARITY. 

C -COMPUTER METHODS FOR MATHEMATICAL COMPUTATIONS-, BY G. E. FORSYTHE, 
C M. A. MALCOLM, AND C. B. MOLER (PRENTICE-HALL, 1977) 

C 

C IN1UT. . 

C 

C NDIM = DECLARED ROW DIMENSION OF ARRAY CONTAINING A 
C N ORDER OF MATRIX 

C A TR I ANGULAR IZED MATRIX OBTAINED FROM SUBROUTINE DECOMP 

C B RIGHT HAND SIDE VECTOR 

C IPVT = PIVOT VECTOR OBTAINED FROM DECOMP 

C 

C OUTPUT.. 

C 

C B = SOLUTION VECTOR, X 
C 

C DO THE FORWARD ELIMINATION. 
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IF (N.EQ.l) GO TO 50 
NM1-N-1 
DO 20 K-l.NMl 
KP1-K+1 
M-IPVT(K) 

T-B (M) 

B (M) -B (K) 

B(K)-T 

DO 10 I-KPl.N 
10 B(I)-B(I)+A(I,K)*T 

20 CONTINUE 

C NOW DO THE BACA SUBSTITUTION. 

DO 40 KB-l.NMl 
KM1-N-KB 
K-KM1+1 

B(K)-B(K) /A(K,K) 

T— B (K) 

DO 30 I-l.KMl 
30 B(I)-B(I)+A(I,K)*T 

40 CONTINUE 
50 B(l)-B(l)/A(l,l) 

RETURN 

END 
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motion is provided by an automatic grinding machine controlled by a computer. The TCA program simulates the 
meshing and bearing contact of the misaligned gears. The transmission errors are also determined. 
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